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Abstrac;t. The small dispersion limit of the focusing nonlinear Scliroodinger equation (NLS) 
exhibits a rich structure of sharply separated regions exhibiting disparate rapid oscillations at 
microscopic scales. The non self-adjoint scattering problem and ill-posed limiting Whitham equa- 
tions associated to focusing NLS make rigorous asymptotic results difficult. Previous studies 
[KMM03, TVZ04, TVZ06] have focused on special classes of analj4ic initial data for which the 
limiting elliptic Whitham equations are well-posed. In this paper wc consider another exa<;tly 
solvable family of initial data, the family of square barriers, ipo(x) = qX[-L,L] for real amplitudes 
q. Using Ricmann-Hilbcrt techniques wc obtain rigorous pointwise asymptotics for the scmiclas- 
sical limit of focusing NLS globally in space and up to an 0(1) maximal time. In particular, we 
show that the discontinuities in our initial data regularize by the immediate generation of genus 
one oscillations emitted into the support of the initial data. To the best of our knowledge, this 
is the first case in which the genus structure of the semiclassical asymptotics for fNLS have been 
calculated for non-analytic initial data. 



1. Introduction and background 

In this paper we study the semi-classical limit of the Cauchy problem for the focusing nonlinear 
Schrodinger equation (NLS) 

(1) ieV't + ^V-xx + V'lV'l^ 
for the family of square barrier initial data 

(2) Mx) = 

for any choice of constants q, L > 0. In this scaling e is a dispersion parameter which we take to be 
small, < e <C 1. Specifically, our goal in studying the semiclassical limit is to derive a uniform 
asymptotic description of the limiting behavior of the e-parameterized family of solutions ^p{x, t: e) 
of the above Cauchy problem as we let e J, for {x, t) in compact sets. 

The focusing NLS equation is of course very well known, describing phenomena in a diverse 
variety of fields, from the theory of deep water waves [BN67] to the study of Langmuir turbulence 
in plasmas [Zak72]. In fact, under suitable assumptions the NLS equation arrises generically in 
he modeling of any weakly nonlinear, nearly monochromatic, dispersive wave propagation. In 
particular, it has been used extensively in nonlinear optics where it has been used to describe the 
evolution of the electric field envelope of picosecond pulses in monomodal fibers [HT73]. In this 
last setting, the small dispersion limit is an appropriate scaling to model the expanding use of 
dispersion-shifted fiber in applications. 
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An important feature of the focusing NLS equation is its modulational instability. The instability 

can be understood in a number of ways, we present here the nonhnear perspective supported 
by Whitham [Whi99]. Without loss of generality we can write the solution of (1) in the form 
ip{x,t) = ^(a;,t)e'^^^'*^/^ with A{x,t) and S{x,t) real. If we assume that A{x,t) and S{x,t) evolve 
slowly with e, then the solution q{x,t) looks locally like a plane wave with amplitude A(x,t), 
wavenumber dxS{x,t), and frequency —dtS{x,t). Writing tp{x,t) in this form and setting p{x,t) = 
\il){x,t)\'^ = A{x,t)'^ and u{x,t) = lmdx^ogip{x,t) = Sx{x,t) we get the following coupled system 
of equations which is completely equivalent to NLS: 

9p _^ dp ^ du ^ 
dt dx ^ dx 
^"^^ du dp du 

dt dx dx 

Assuming that the wave parameters evolve slowly with e we can regard the right-hand side as 
a perturbativc correction which we may formally drop. Doing so we arrive at the (genus zero) 
modulation equations associated with (1), 
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This is a quasilinear system of partial differential equations and it is easy to check that its char- 
acteristic speeds are given by u ± i^/p which are complex at any point at which p ^ 0. Thus 
the limiting system of equations is elliptic for which the Cauchy problem is ill-posed. In fact, the 
Cauchy problem for an elliptic systems is solvable (by the Cauchy-Kovalevskaya method) only if we 
insist that the initial data p{x,0) and u{x,Q) are analytic functions of x. This is an overly restric- 
tive condition for a physical system, and one which our square barrier initial data (2), with jump 
discontinuities, does not satisfy. The ill-posedness of the limiting equations makes the semiclassical 
limit a very delicate question in general and particularly delicate in the case of interest where the 
initial data is discontinuous. In fact, our results imply that even away from the discontinuities 
the limiting behavior cannot be described by a guess of the form ij) = Ae^^/'^ without admitting 
multi-scale expansions of A and S from the onset. 

The original observation was that the modulational instability of the NLS equation would lead 
to ostensibly chaotic wave-forms. However, a number of numerical studies [MK98, BK99, LM07] 
observed that the wave forms show a large degree of order exhibiting disparate regions of regular 
quasiperiodic oscillatory e-scaled microstructures separated be increasingly sharp curves in the 
(a:, t)-plane called caustics or breaking curves, whose scale is set by the initial data. Our main goal 
in studying the semiclassical limit is to analytically find the caustics and to describe the asymptotic 
behavior of the oscillations in the interlaying regions offset by these caustics. 

Analytic results concerning the behavior of NLS are possible because of its integrable structure 
which allows one, in principle, to solve the Cauchy problem by the forward/inverse scattering 
transform using the Lax pair representation [ZS71] for NLS. This procedure consists of three steps: 

Step 1. Prom the eigenvalue problem in the Lax pair (5a) determine the "spectral data" produced 

by the given initial data t/jq. This is the forward scattering step in the procedure. 
Step 2. Use the time- flow part of the Lax pair (5b) to describe the spectral data's time evolution 

as the potential V' evolves under the NLS equation. 
Step 3. Given the scattering data at time i, invert the forward scattering transform of Step 1 to 

find the solution ip{x, t) of the NLS equation at a later time. This step for NLS, as is often 
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the case, may be characterized as a matrix-valued Riemann-Hilbert problem, reducing the 
inversion problem to a question of analytic function theory. 

For a non-integrablc problem there is no reason to believe such a procedure should exist. The 
amazing fact is that for integrable problems the evolution of the spectral data in step 2 is explicit 
and relatively simple, and moreover, the inversion of the spectral map required in step 3 is actually 
a well-defined map. 

The semiclassical limit greatly complicates the analysis in both Steps 1 and 3 In the scattering 
step one must determine both the discrete spectrum (i^ eigenvalues) and the so called reflection 
coefficient defined on the real line (the continuous spcctriim) of the Zakharov-Shabat (hereafter 
ZS) scattering problem (5a). The focusing ZS problem, unlike the scattering problems for KdV or 
defocusing NLS, is not self-adjoint and this makes the determination of the scattering data highly 
nontrivial. Numerical evidence [Bro96] and asymptotic results [MilOl] suggest that the number of 
eigenvalues scales as 1/e and asymptotically accumulate on a "Y" -shaped spectral curve for a large 
class of initial data. The analysis of the inverse scattering step is also frustrated by the vast amount 
of spectral data. Among other issues, the eigenvalues of the scattering step become poles of the 
Riemann-Hilbert problem and one is faced with seeking a function with an ever growing number 
of poles. The question of how to best interpolate the poles is one of the first difficulties one has to 
address in the inverse scattering step. 

The situation is somewhat simplified for real data with a single maximum, in this case it has been 
shown [KS03] that the eigenvalues are confined to the imaginary axis. However, exact determination 
of the eigenvalues is possible only in special cases, among them the Satsuma-Yajima data and its 
generalization [SY74, TVOO] , and piecewise constant data. The exact solvability of the scattering 
problem for square barrier data (2) is one of the principle motivations for using it as a model 
of how the semiclassical evolution smooths out discontinuities. The somewhat common practice 
of replacing the exact scattering data with simpler approximate data is a dubious procedure for 
studying initial data with discontinuities. The scattering problem has been shown to be very 
sensitive to perturbations of the initial data [BroOl], and its not clear how using "nearby" scattering 
data would effect the discontinuities of the initial data. For this reason we found it crucial to solve 
the true Cauchy problem; we make no approximations to either the initial or scattering data. The 
only approximations we make are in the inverse problem where the error can be controlled. 

Despite the numerous frustrations, in recent years there has been considerable progress on study- 
ing the semiclassical focusing NLS equation. In [KMM03] the authors considered the semiclassical 
limit of certain soliton ensembles, that is, families of real, analytic, reflectionlcss initial data whose 
eigenvalues are given by a quantization rule. In a separate study [TVZ04, TVZ06] the initial data 
ipoix) = sech(a;)^"'"*^/^ was considered to study the effect of a nontrivial complex phase on the evo- 
lution. In both cases, the authors found sharply separated regions of the (.T,t)-plane inside which 
the leading order asymptotic behavior of the solution is described in terms of slowly modulated 
wave trains of a particular genus whose Riemann invariants' evolution is described by the Whitham 
equations of the corresponding genus. For each fixed x and increasing t one encounters these sharp 
transitions, or caustics, across which the genus of the asymptotic description changes. In particu- 
lar, in both cases it was observed that up to some primary caustic T{x), that is for t < T{x), the 
evolution was initial described by slowly modulated plane wave with invariants satisfying the genus 
zero Whitham equations (4). 

In both the above cases the initial data considered is analytic, and this plays an important role in 
the analysis of the inverse problem. Additionally, we note that due to the analyticity of the initial 
data the limiting elliptic modulation equations are locally well-posed. For non-analytic initial data 
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it's not clear what should happen for small (but fixed) times in the semiclassical limit; in principle, 
the evolution could be described by arbitrarily high genus waves from the onset of evolution. The 
situation for square barrier data is necessarily complicated; in the small time limit it was shown first 
in [DM05] for the defocusing case and in [ER06] for the Manakov system (of which focusing NLS 
is a special case) that the small time limit of the evolution with initial data given by (2) is equal 
to the evolution of the linear problem and Gibbs phenomena occur at the discontinuities. These 
Gibbs phenomena generate oscillations of arbitrarily short wavelength which propagate away from 
the discontinuities. Of course, in the semi-classical scaling any fixed time independent of e is not a 
small time in the limit, but the existence of rapid oscillations from the onset is a strong indication 
that something singular must be occurring near the discontinuities. One of the principal goals in 
this study is to understand how the semiclassical NLS evolution smooths the discontinuities in the 
square barrier (2) initial data. 

In this paper as in all of the above mentioned studies the bulk of the work is contained in 
analyzing the Riemann-Hilbert problem (hereafter RHP) associated with the inverse scattering 
transform. The results are asymptotic in nature and rely on the nonlinear steepest-descent method 
developed by Deift and Zhou first in [DZ93] and then generalized in [DZ94, DVZ94] and in the 
works of many others. The method, in a nutshell, seeks to construct an explicit parametrix P{z) 
such that the transformation m{z) = E{z)P{z) maps the solution m{z) of the original RHP to 
a new unknown E{z) which satisfies a small-norm RHP whose solution and asymptotic behavior 
can be given in terms of certain singular integrals. In each of the semiclassical and small-time 
studies mentioned above, as in most cases where the steepest-descent method has been successful, 
the construction of the parametrix revolves around the study of a single complex phase function of 
the type e*^^^'/'' which must be controlled. 

The fact that the scattering data in the previous studies can be reduced to a single complex phase 
is a consequence of the choice to study very 'nice' initial data. By studying square barrier data we 
hope to gain an understanding of solutions for more generic types of initial data, i.e. data outside 
the analytic class for which the limiting problem (4) is ill-posed. As soon as we consider initial 
data outside the analytic class, we immediately encounter scattering data which cannot be reduced 
to a single dominant complex phase. For the square barrier data considered here the reflection 
coefficient takes the form of a whole series: X^^^q where the rk{z) are e-independent 

algebraically decaying weights (c.f. Sec. 3.1.1). As we will show, as x and t evolve, different 
terms in the expansion give the dominant contribution in different parts of the spectral z-plane. 
In the course of our analysis we show that for times t up to the secondary breaking the analysis 
is dominated by the two lowest phases 6q and 6\. However, we see no a-priori way to rule out the 
possibility that at later times the higher order harmonics emerge to contribute at leading order. 

In addition to being multi-phased, the reflection coefficient for rough initial data is necessarily 
slowly decaying; just as the Fourier transform the reflection coefficient, r{z), for discontinuous initial 
data decays as 1/z for z G K. The slow decay forces us to consider the contribution of the reflection 
coefficient everywhere on the real line. This makes the analysis more difficult than those considered 
in [KMM03, TVZ04] where the reflection coefficient, at worst, contributes only on a finite interval. 
To complete the inverse analysis we had to marry techniques for analyzing multi-soliton RHPs with 
those for pure reflection problems. This complication is greatest at the real stationary phase points 
introduced by the steepest-descent method. On the real line no phase is completely dominant and 
the multiphase nature of the reflection coefficient lead us to construct delicate 'annular' local models 
in neighborhoods of these stationary phase points. This model is discussed in Section 3.3.3; to our 
knowledge our annular construction is the first of its kind. 
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1.1. Statement of the main results. The culmination of our work in this paper is the foUowing 
theorem which describes the leading order asymptotic behavior of the solution of the Cauchy prob- 
lem defined by eqs. (1) and (2) for arbitrarily large times t outside the initial support and up to a 
secondary caustic for x inside the support of the initial data. 




Figure 1. The regions Sk, k = 0,1,2 in which we have a description of the 
semi-classical limiting behavior for the initial data (2) and the caustics Ti and T2 
separating the regions. 



Theorem 1. Consider the curves t — Ti(|a;|) andt = T2(|a;|) defined by (89) and (144) respectively. 
Label by Sk, k = 0,1,2, the following subsets of x [0,cx3) (see Figure 1): 



Sq = {{x,t) 
Si^{ix,t) 
S2 = {{x,t) 



\x\ > L,t>0} 
\x\<L,0<t< Ti{x)} 
\x\<L,T,{x)<T<T2{x)} 



Then for all sufficiently small e, the solution rpix, t) of the focusing NLS equation (1) with initial date 
given by (2) satisfies the the following asymptotic description on any compact Kj d Sj, j = 0,1,2: 



mx,t)-i;asyix,t)\^o[e'^Hoge^ , 



where, 
i. For {x, t) G So, 



1pasy{x,t) = 0. 



ii. For {x,t) £ Si, the limiting behavior is described by the nearly plane wave solution 



'4'asy{x,t) = gexp 



{q'^t + eu!{x,t)) 



where 



uj{x, t) = — 



1 



log(l-K|ro(A)|2) 



2 ^ql 



dA, 



andro, ^i, and are defined by (29), (90), and (92) respectively. 
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Hi. For {x,t) e S2, the limiting behavior is described by the slowly modulating one-phase wave-train 

ipasy{x,t) = {q- Ima))- 



9(0) B(2_l(^) + /To-K-ifi) 



e(2A(oo)) e(iTo - ie-^il) 

Here, 6(z) is the Reimann theta function, see (170), for the Riemann surface associated with 

the function 

H{z) = \/{z - iq){z + iq){z - a{x,t)){z - a{x,t)*) 

The points a{x, t) and a{x, t)* are complex Riemann invariants described by Propostition 9 and 
whose evolution satisfies the genus-one Whitham equation (141). The remaining parameters 
are all slowly varying functions of {x, t) described in terms of certain Abelian integrals: 



, . r &z / r dz 



and O, To, and Yq are given by (125), (174), and (177) respectively. 

We make the following observations concerning our results: 

1. The barrier does not asymptotically spread out in time, though at any time t > the solution 
is no longer compact supported. Our theorem establishes an asymptotically vanishing bound on 
the amplitude of the solution for x outside the support of the initial data. However, it does not 
provide any information about the nature of the small amplitude oscillations that occur in the 
exterior region. These smaller corrections could be found by calculating the higher order terms 
in the expansion of the small norm RHP 3.5 but we did not pursue this here. 

2. In the regions So (trivially) and iSi the asymptotic limiting solution tpasy{x,t) is a slowly mod- 
ulating plane wave solutions of (4), the genus-zero Whitham equations. The most interesting 
feature of the solution in this region is the slowly evolving correction to the phase uj{x,t). This 
correction has in important consequence: the Whitham modulation theory (discussed around 
(3) and (4) above) fails to capture the behavior of the solution. The modulation theory cannot 
recover the slow phase correction Lo{x,t) which constitutes an O (1) correction to the solution's 
phase; the Riemann-Hilbert analysis naturally recovers this phase. This is explained in detail in 
Section 4.7 where we show that this correction does not correspond to a regular correction to 
the geometric optics approximation predicted by WKB. 

3. The evolution regularizes the initial shocks in the square barrier data at x = ±L by the instan- 
taneous onset of genus-one oscillations in the growing region 1S2. This is consistent with the 
previous studies [KMM03, TVZ04] where the mechanism for formation of the primary caustic is 
the formation of an elliptical shock in the solution of the genus zero modulation equations valid 
for times before the primary caustic. 

4. The fact that the reflection coefhcient is NOT single phase is not merely a technical difficulty. 
The second breaking time T2 (x) which gives an upper bound on the genus-one region S2 can be 
characterized in terms of a critical transition related to the second phase. 

1.2. Outline of the rest of the paper. The remainder of the paper is concerned with establishing 
the proof of Theorem 1. In Section 2 we quickly review the details of the forward scattering theory 

pertinent to our analysis and show how one explicitly calculates the scattering data for initial 
data given by (2). In section 3 we consider the inverse problem for x outside the support of the 
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initial data. Section 4 concerns the inverse problem inside the support up the first caustic while 
Section 5 covers the inverse problem between the first and secondary caustics. The analysis in 
each of the last three sections is largely self-contained though for the sake of brevity we use results 
and estimates from the previous sections whenever possible. In particular, the properties of the 
reflection coefficient and resulting jump matrix factorizations discussed in Sections 3.1.1 and 3.2 
are used throughout the paper. 

Finally, a quick note on notation. Throughout the paper we make extensive use of the spin 
matrices: 

^0 1\ /O -i\ [I 



Additionally, given any nonzero scalar A and matrix A let 



Q , and A*^'^M = A'^MA-'^^ 

With respect to complex variable notation, we let z* denote complex conjugate of any complex 
number or matrix z and let A'^ denote the hermitian conjugate of any matrix A. For complex 
valued functions f{z) of a complex variable z we denote by /* (z) or simply /* the Schwarz reflection: 
f{z*)*. In the rare case when only the function value is to be conjugated we write f{z)*. 

2. Forward scattering of square barrier data 

We here quickly review the forward scattering theory necessary to define the Riemann-Hilbert 
problem associated with the inverse scattering problem. We omit proofs, preferring to get as quickly 
as possible to the analysis of the inverse problem. The interested reader is referred to [BC84, BDT88] 
for a comprehensive treatment of the subject. 

The NLS equation is completely integrable in the sense that it has a Lax pair structure. In 
[ZS71] Zakharov and Shabat showed that the focusing NLS equation (1) can be recognized as the 
compatibility condition for the following overdetermined pair of equations: 





(5a) eW^ = jCW, C := -iza^ 



-^p* 



(5b) ieW. = BW, B:=izL-\{^l 

The Lax pair allows one, in principal, to solve the NLS equation using the forward/inverse scattering 
procedure. In the forward scattering step, one takes given initial data 'ijj{x,t = 0) = tpoi^^) that 
decays suitable quickly and seeks a solution W of the first equation in the Lax pair (5a), known as 
the Zakharov-Shabat eigenvalue problem, in the form 

(6) W = me-i='"''' 

such that, 

i. m ^ I, as a; ^ 00 

(7) 

a. rn {-, z) is bounded and absolutely continuous 

The construction of m from initial data is a nontrivial step which can be calculated explicitly only 
in special cases. However, for any potential Vo € ^^(M), the solution m must have the following 
properties [BC84]: 

1. There exist a discrete set Z in C\M such that for each z € C\(RUZ) there exist a unique solution 
of W of eW^ = CW in the form (6) satisfying (7), 
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2. m{x, z) is a meromorphic function of z in C\M whose poles are precisely the elements of Z and 

on C\IR, lim^^oo 'rn{x, z) = /. 

3. The function m{x,z) takes continuous boundary values m±{x,z) := lim5_).o+ m{x,z ± id) for 
each 2; e K which satisfy a "jump relation", m+ = m-v{x,z), where: 

(8) + 1^(^)1' 



r{z) 1 

4. The points in Z come in complex conjugate pairs. At each simple pole (the generic case) the 
residues of m satisfy a relation of the form 



for Zk € C"*", Resm = limm(a;, z) 

(9) 

for zl G C~, Resm = limm(a;, z) 





Cke~^ 



-c*e-^ 




The form of the jump matrix and the poles coming in conjugate pairs is a consequence of the 
reflection symmetry 

(10) m{x,z*) = a2m{x,z)* (72- 
implied by (5a). 

The collection of poles {zk}^ the connection coefficients {c^}, and the function r{z), called the 
reflection coefficient, appearing in the jump matrix (8) constitute the totality of the scattering data 
generated by the initial potential tpoix). 

Though one is interested in building m(x, z) for z off the real axis, the construction of the 
scattering data usually bc;gins by restricting z to the real axis. For real z one seeks solutions W^~^^ 
and of (5a) normalization such that: 

(11) lim W(±)ei^^'^^ = I. 

X— )-±oo 

Each of these so called Jost solution constitutes a fundamental solution of the differential equation 
for 2: G M and thus the two solutions satisfy a linear relation of the form 

(12) W'^-\x,z) = W'^+\x,z)S{z), zeR. 

The matrix S, called the scattering matrix, is independent of x and takes the form 

(13) S{z) = (jg , det 5 = \a{z)\' + \b{z)\' ^ 1. 

The normalized solution m is then constructed by examining the integral equation representations 
of W^^^ and observing that different columns of each solution extend analytically to C+ or C^. In 
the process of this construction one finds that the function a{z) extends analytic to C+, while b{z) 
is, in general, defined only on the real axis. Moreover, the poles Zk are precisely the zeros of a{z) 
in C+ and the reflection coefficient r{z) is defined by the ratio 

(14) r{z) :=b{z)/a{z), zeR. 

The time evolution of the scattering data is remarkably simple. As ijj evolves according to (1) 

the coefficients in the Lax operators C and B gain time dependence which in turn is inherited by the 
spectral data. A basic consequence of the Lax pair structure, is that the time-flow is iso-spectral 
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[Lax68], that is, the eigenvalues {zk} and their conjugates are invariant in time. The remaining 
scattering data evolve as follows: 



2itz 



(15) r{z,t)=r{z)e « , and Ck{t) = Cke' 

where r{z) and C). above arc the values of each datum at time t = 0. 

Given the time evolved spectral data, the problem of reconstructing the potential ijj can be 
cast as a meromorphic Riemann-Hilbert problem (RHP). For each x and t one seeks a piecewise 
meromorphic matrix-valued solution m(z; .x, t) of the following properties: 

Riemann-Hilbert Problem 2.1 for focusing NLS. Given the intitial scattering data r{z), 
{zk}, and {ck}, find a function m{z;x,t) such that: 

1. m is meromorphic in C\M with poles at the points Zk and their complex conjugates. 

2. As z ^ oo in C\K, m I . 

3. As z approaches K from above and below m assumes continuous boundary values m± satisfying 
the jump relation m+ = ru-V, where 

(16) V = e-i(*^^+-) ei(*^^+-). 

^. The poles of m at each Zk and z^ are simple and satisfy the residue conditions 

0^ 



^^'^'^ /q _^*^-'-{2tz^+2xz) 

Res m = lim " 





If m{z;x,t) is the solution of this problem for the giving scattering data generated by initial data 
■00 (x) then the function defined by the limit 

(18) i^{x,t) = 2i lim zmi2{z;x,t) 
exists and is a solution of the focusing NLS equation (1). 

2.0.1. Scattering for compactly support data, exact solutions for square barrier data. Here we record 
some additional properties of the scattering data when the initial potential ipo{x) is compactly 
supported. Additionally, we find explicit formulas for the scattering data when the initial data is 
piecewise constant and compactly supported which includes the square barrier initial data with 
which we are interested. 

For bounded compact potentials ipQ, the Jost functions W'^^''{x, z) are entire functions of z. The 
scattering relation (12) thus holds not only on the real axis, but in the entire complex plane, and the 
scattering coefficients a{z) and b{z) are both entire functions for compact initial data. Additionally, 
the connection coefBcients Cfe can be expressed explicitly in terms of the scattering coefficients as 

(19) c, = -^=Resr(^). 

a'{zk) 

This relation between the connection coefBcients and the residues of the reflection coefBcient, which 
is false for generic initial data, will later play a fundamental role in our analysis of the inverse 
problem. 

As discussed previously, calculating explicit formulae for the scattering data generated by generic 

initial data is intractable, particularly in light of the singular dependence of the scattering data on 
the dispersion parameter e. However, for initial data which is compact and piecewise constant the 



10 



ROBERT JENKINS AND KENNETH D. T-R MCLAUGHLIN 



forward scattering problem is simple. We can assume freely that the initial data is supported on 
an interval [—L,L], and write it as 

n 
k=l 

where —L = Xq < Xi < ... <a;„ = Z/isa partition of the interval [— L,L], the qk are complex 
constants, and \a,b) denotes the characteristic function on [a, 6). On each interval the 

ZS scattering problem (5a) reduces to a constant coefficient ODE. The Jost solutions W'^^^x^ z) 
are calculated by solving the ZS problem on each interval and patching the solution together by 
demanding continuity at the partition points. That calculation leads to the following formula for 
the scattering matrix: 

Q\{-izcr3+Qn){Xr,-X^-l) ^izL(T3/€ Q ( ^ 9k 



-<il 0^ 

In the simplest case, the potential consist of a single barrier Vo(2;) = ?1[-l,z,]- The scattering 
coefficients a{z) and h{z) are then 

viz) COS — — iz sm — 

(20) a{z) =— ^-^^^ ^ ) ^2Liz/e^ 

2Lv(z)'' 



—q sm 
'^'^ = —.iz) ' 
where 

(22) ly = + q2_ 

From the explicit formalae it is clear that both a(z) and b{z) are both entire functions of z. For 
real values of the constant q, the initial potential ^0(2;) falls in to the class of real potentials of a 
single maximum; as such the poles Zk, which are the zeros of a{z) in the upper half plane, must lie 
on the imaginary axis [KS02]. Using the explicit formula for a{z) wc find that the zj. are confined 
to the imaginary interval i[0,q), and asymptotically there are roughly 1/e poles in the the interval; 
as e 4, they collect according to the asymptotic density 

2Lq z 



(23) po{z) = 



which shows that the density has intcgrablc singularities at z = ztiq. For convenience we exclude 
the countable set of values at which eigenvalues are 'born' at the origin 

_ 4Lq 

^ r ^ni — 77; ; ) — 0, 1, 2, ... . 

(2n + IjTT 

Making use of these explicit formulae for the scattering data wc state the exact Riemann-Hilbert 
problem for the inverse scattering given square barrier initial data defined by (2); 

Riemann-Hilbert Problem 2.2 for focusing NLS with square barrier initicil data Find 
a matrix valued function m{z; x, t) satisfying the following properties: 

1. m is meromorphic for z € C\K, with simple poles at each Zk € C"*" and its complex conjugate. 
Here, the points Zj~ enumerate the zeros in C"*" of the function a{z) defined by (20). 
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2. ^5 2; oo, m = 7 + 0(1/2). 

3. m takes continuous boundary values for z G 
rri-V where 



I, m+ and rri- satisfying the jump relation, m+ 



(24) 



V = 



re 



4- m satisfies the residue conditions 



6 = 2tz^ + 2xz. 



(25) 



Res m = lim m 



Res r*e 




In (24) and (25), the reflection coefficient r{z) is given by 



(26) 



r(z) 



-gsm 



v\z) COS ( — ^ \ — iz sin ( — ^ 1 



-2Liz/t 



Finally, if m{z;x,t) is the solution of this problem, then (18) gives the solution of focusing NLS 

(1) given the square barrier initial data (2). 

The goal now is to use the Deift-Zhou steepest-descent method to construct an approximation 
P{z) to the solution of this RHP such that the ratio E{z) defined by m{z) = E{z)P{z) is uniformly 
accurate in the scmiclassical limit as e — > 0^. The approximation depends parametrically on [x^t) 
via the function 9{x,t) appearing in the RHP. In the following sections we will construct such an 
approximation in each of the space time regions Sk, k = 0,1,2, described in Theorem 1. This 
procedure consists of four major steps. 1. Wc introduce a pole removing transformation m h- > Af 
such that the new unknown M satisfies a holomorphic RHP. 2. The transformation to M typically 
introduces contours on which the jump matrices are asymptotically unstable in the scmiclassical 
limit; wc therefore introduce a so-called g- function transformation to remove the unstable jump 
matrices. 3. By explicit factorization of the jump matrices we deform the remaining oscillatory 
jumps onto contours where they are exponentially near identity. 4. We build a global approximation 
A{z) to the remaining problem and show that the resulting error E{z) satisfies a small norm RHP. 
Once the problem is reduced to small-norm form, the theory of small norm RHPs gives a uniform 
asymptotic expansion of E{z). By unravelling the series of explicit transformations leading from 
TO to wc get an asymptotic expansion for the solution of the original RHP 2.2 and through (18) 
the solution ^{x, t) of NLS. 



3. Inverse problem outside the support of the initial data 

Here we consider the Riemann-Hilbert problem 2.2 with scattering data generated by the square 

barrier initial data (2) for values of x outside the support of the initial data. Because the NLS 
equation preserves even initial data we will consider only x > L. As will become apparent this 
is the simplest case in which to carry out the inverse analysis and much of what we do here will 
be used in Sections 4 and 5 where we consider the inverse analysis for x inside the support of the 
initial data. 
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3.1. Reduction to a holomorphic RHP. To utilize the nonlinear steepest descent method we 

first need to remove the poles from the Riemann-Hilbcrt problem. In the rcflcctionlcss case this can 
be accomplished by introducing contours C and C* enclosing the locus of pole accumulation in C+ 
and C~ respectively. A new unknown M is constructed inside these contours from an interpolate 
of the connection coefficients and an explicit Blaschkc factor term encoding the Zk in such a way 
that the new unknown has no poles, see [KMM03] for such a construction. For compact initial data 
there is a more direct construction; the reflection coefficient extends meromorphically off the real 
line and the relation (19) shows that the reflection coefficient naturally interpolates the poles of m. 



1.5 

1.0 

05 

0.0 
-0.5 
-1.0 
-1.5 

-1.0 -0.5 0.0 05 1.0 1.5 2.0 

FiGUR.E 2. Schematic representation of the; contours C and C* and the regions f7 
and n* involved in the pole removing factorization. The dashed line represents the 
loci of pole accumulation. 




Fix a point ^ < 0, and take C a semi-infinite contour in C"*" terminating at ^ such that the region 
n enclosed by C and the real interval [^,oo) contains the locus of pole accumulation 7'[0, g], see 
Figure 2. Let C* and il* denote the respective complex conjugate contour and region. Then the 
new unknown 



(27) 



H 



1 

1 r*e- 



m 



: z eSl 

: zeQ* 
: elsewhere 



has no poles. The price one pays to remove the poles in this way is the appearance of new jumps 
along the contours C and C* in the resulting RHP for H. In order to preserve the normalization 
of the RHP and arrive at a problem amenable to steepest descent analysis the transformation 
from m to H needs to be asymptotically near-identity for large z. This condition amounts to 
understanding how the factor re'^/* appearing in the off-diagonal entries of the transformation 
behaves in the complex plane. 
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3.1.1. Multi-harmonic expansion of the reflection coefficient. For z £ M the reflection coefficient 
(26) is rapidly oscillatory, but once z moves off the real axis any finite distance it has a simple, 
e independent, asymptotic limit. As previously noted, the values of the reflection coefficient are 
independent of the branch cut given to v ~ + q'^. For convenience we take v to be branched 
on some simple finite contour 7 and normalized to ^ z for z large. For z G C+ outside the region 
enclosed by 7 and the imaginary axis (i.e. the set where Im v > 0) 

00 

(28) r(z)e^«/^ = ^rfe(z)e''''=/% 

fc=0 

where 
and 

rfc(z) = (-l)^-Vo(z)2'=-i(l + |ro(z)p), fc > 1; 

9k^2tz'^ + 2{x- L)z + 4:kLiy, k>0. 

In particular, the expansion is uniformly convergent for z e with \z\ large. The expansion 
highlights the central technical difficulty of this problem; in previous studies where the semi-classical 
inverse problem has been successfully solved [KMM03] , [TVZ06] , a single exponential factor emerges 
in the jump matrices which one must control. Here, we have an infinite sum of different harmonics — 
one can think of (28) as a generalized Fourier series for r — each contributing to the analysis. A 
priori, we have no way of knowing which harmonics contribute for each choice of parameters {x,t). 
However, in the course of our analysis we will see that, at least up to the second caustic, the first 
two harmonics 9q and 9i dominate the analysis. 




Figure 3. The regions of growth (white) and decay (grey) of the two principal 
harmonics for x > L and t > Q. The dashed line represents the locus of pole 
accumulation for m{z). 



The pole removing factorization (27) works if the region fi, except for possibly a compact subset, 
lies in a region in which exponentially decaying. This amounts to understanding for each 

Ok where Im^^ is positive (decay) or negative (growth). For a; > L we fix 7, the finite branch cut of 
V, to lie on the imaginary axis. This choice of branch makes Im v{z) > for all z in C+\ «(0, q\ and 
thus \Ta.Ok{z) is an increasing function of k for each fixed z; if we let D}. = {z £ C+ : Imflfc > 0} 
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then -Dfc-i C -Dfe for A; > 0. Furthermore, it is easily shown that each 6k has a single, real-valued, 
stationary phase point ^fe, defined as the unique solution of j^Ok = 0, and as k increases, the ^^s 
decreases. 

The regions of exponential decay and growth of e'^°/^ and e'^^/* for generic values of a; > L and 
t > are plotted in Figure 3. The depend parametrically on x and t and we have 

(31) f«o) = . = 

For fc > 1 no such simple formula for the can be given, however, for small times we have 

x+{2k-l)L 



(32) 



2t 



+ 0{t). 



In particular, for x > L, < and the locus of pole accumulation in the upper half plane, i[0,<7], 
lies entirely in the region of exponential decay of all the 0k- Motivated by these observations, let 
Fo be a semi-infinite contour in C+ emerging from such that it encloses the interval i[0, q], stays 
completely within the region of exponential decay of e*^°/^ and is oriented away from ^o; let r^o 
be the region between Fq and the positive real axis; let Fg and Qq be their respective complex 
conjugates, see Figure 5. Define 



(33) 



1 

— re«' 



M = ( 



m 



re 



z e 

elsewhere 



As observed previously, it follows from (25) that the new unknown M will have no poles. It follows 
from its definition that M satisfies the following problem. 

Riemann-Hilbert Problem 3.1 for M: Find a matrix valued function M such that: 

1. M is holomorphic for z € C\Fm where Tm = Fq U Fq U (— oo, ^q]- 

2. As z oo, M = 1 + {1/z). 

3. M takes continuous boundary values for z G Fm, M+ and M_ satisfying the jump relation, 
M+ = M_Vm where 



(34) 



Vm = < 



1 



re 
1 



j,*g-if/e 
1 



1 



z e (-oo,Co) 



z€To 



The new jumps follow from calculating (M_)~^M+ on each contour. Observe that the new 
jumps on Fq and FJ are exponentially near identity away from the real axis. Thus for x > _L, the 
poles can be removed without introducing any poorly conditioned jumps; as we will see in Sections 
4 and 5 this is markedly different from the situation for x G [0, L) and accounts, to large extent, for 
the difference in the resulting asymptotic description of the solutions inside and outside the initial 
support. 
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3.2. Reduction to Model Problem, deformation of the jump contours . Now that the 

poles have been removed, wc are left with a RHP with rapidly oscillatory jumps on (— 00,^0] • To 
employ the steepest descent method we seek by explicit transformations to deform these jumps to 
regions where they are exponentially decaying. Our principal tool for finding such transformations 
is matrix factorization. Given an oscillatory jump matrix f on a contour C, we seek factorizations 
of the form v = b-vb'^^ such that the factors b± are analytically extendable off C. If we then 
introduce contours C+ and C_ , to the left and right of C respectively, see Figure 4, we can define 
a new unknown 



rrir, 




z e Int(CUC_) 
z e Int (C U C+) 
elsewhere 



The new unknown m„eu) acquires new jumps Vnew = {ninew)-^{'mnew)+ equal to 




Figure 4. Given a factorization v = b- vb'^ of the original jump matrix on C, 
the definition of a new unknown (left) and the resulting new jump matrices (right) 
in a generic lens opening. 







zeC 


. -J 

'-'new — 




zeC+ 






zgC_ 



Such a factorization is useful provided that and &_ are near identity on their respective contours 
and the jump remaining on C, v, is no longer rapidly oscillatory. 

3.2.1. Factorization to the left 0/^1. The original jump matrix v (cf. (24)) has two common 
factorizations: 

(35) V = R^'R, 

(36) V = R{l+rr*f^ R''. 
Where 



(37) R:={^^ig/, J and R-=y_^^i0/e ^ 



The first factorization was used already in the factorization that removed the poles from the original 

RHP. For Re z < ^1 we are interested in the second of these factorizations. The rightmost factor 
W is the factor which will factor into C"*", we need to understand the behavior of its off-diagonal 
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entry in the upper half-plane. Recalling the definition (14) of the reflection coefficient and the 
unimodularity (12) of the scattering matrix we have 



1 + rr* = 1 + 



bb* 



1 

aa* 



= b*a. 



1 + rr* 

Using (20) and (21) we can express a and b* in the form 



(38) 



6* 



z + 1/ 

iq 



,2Li{z-v)/t 



-2Liv/t 



1 



2 ALiu/e 



and these expressions give the off-diagonal entry in W in the form 



(39) 



1 + rr* 



-ie/t _ 



(1 



1 _ giLiv/e 



1 - rle^"^'"/' 



Recalling that lnw > for each z G C+\2(0,q] it follows that W is exponentially near identity in 
any region in which the exponential e*^i/'= is growing. Thus, (36) is a good candidate factorization 
for Rez < ^i. The only issue is that the central diagonal factor (1 + rr 



to be left on the real 



axis is still rapidly oscillatory as e 4 0: this term must also be factored. Using (38) define 



(40) 



ao := 



z + v 



•2Li{z-v)lt _ 



J2Li{z-v)/e 



which captures the leading order behavior of the function a for z G C+. Note that for z G C"'"\(M U 
(0, ig]) Imi/ > so the quantity a/oo = 1 — rge^^"'/^ is exponentially near one. Motivated by this 
observation we introduce the factorization 



(41) 



v={R{a*/al)-''-) . (1 + rorS)^^^ . ((a/ao)-^=^^t) 



factors into C~ 



stays on ] 



factors into C+ 



of the original jump (24) on the interval (— oo, ^i). This factorization greatly simplifies the analysis 
of the RHP; it reduces the problem to a consideration of only the first two harmonics and 9\ in 
the expansion (28) of the reflection coefficient. 



3.2.2. Factorization on the interval ($i,^o) • This interval is the boundary on the real axis of an 
open region of C"*" in which only the zero*'' order harmonic e'^"/^ is large, see Figure 3. Motivated 
by (28) and the matrix factorization (35) we introduce the factorization 



(42) 



roe 



■ roTo roe 

ieo/e 



By Rq, and later _Ro, we denote matrices of the same form as R and R respectively, where in the 
off diagonal entries re*^/^ is replaced by its leading order approximation 



(43) 



roe 



1 



and i?o 



1 



l+ror-j; 
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The exterior factors RRq ^ and its hermitian conjugate in (42) 'subtract out' the dependence on 



Q^Oo/e^ in the (2, l)-entry we have 

(44) re 



ie/e , 



roe 



ieo/e 



ro 



1 - r^e'^^^'^A 



1 - rle'^^^"/'' 



which at leading ordcT depends on e^^^ ^'^ which is asymptotically small in for Re z > . This 
allows us to open lenses (see (46)) which deform these exterior factors into the regions il\ and 
depicted in Figure 5. 

The remaining middle factor in (42) takes the same form as the original jump matrix with the 
reflection coefficient replaced with its leading order behavior. In the region of the upper half-plane 
bounded by and the harmonic e'^"/^ is exponentially large. This factor is deformed off the 
axis by an analog of three term factorization (36) (the same terms with re^^l^ replaced by roe'^°/*) 
yielding the final factorization 



(45) 



V = [r^Ro^Ro) . (1 + ror^or ■ (rIRo'R) • 



factors into C" 



stays on ] 



factors into C+ 



Fi 




Figure 5. The contours and regions flk used to define the lens opening trans- 
formation M Q. Their exact shapes are not important, each is chosen to lie 
in regions where the corresponding factorizations given in (46) are asymptotically 
near identity. Note, that the contour Fa is oriented from to ^o- 

3.2.3. Reduction from a holom,orphic RHP to a stable model problem. Using the factorizations given 
above we are ready to define a new unknown Q{z). To define Q we introduce the following set of 
contours and regions and their complex conjugates. Let Fi and F2 be rays in C"*" meeting the real 
axis at ^1 such that away from the axis they remain completely in regions where e^^i/*^ is decaying 
and growing respectively. Additionally we take Fi to be bounded away from Fq. Let F3 be a contour 
in C+ terminating at ^0 and ^1 such that it lies entirely in the region where e'^°/^ is growing and 
does not intersect either Fq or Fi. Finally, let fli be the region bounded between Fq, Fi, and F3; 
similarly let ^2 and fl^ be the region between the real axis and F2 or F3 respectively. The contours 
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Ffe, their orientations, and the sets Clk are depicted in Figure 5. Define 



(46) 



Q 



'MRoR-^ 




MR-^{a/aoY^ 


Z&Q.2 


MR-^RoRq'' 


z&Q.^ 


MR^Rq^Ro 


z 


MR{a*/a*o)-''3 


Z &Q.*2 


MR^Rq^ 


Z&Q.I 




elsewhere 



Let Fq denote the union of contours used to open lenses, Fq = Ufc=o(-'^fe ^X) denote 
the totality of contours on which Q is non-analytic, Tq = Tq U (— oo, ^o] • Using the factorizations 

(41), (42), and (45) it follows that Q satisfies 

Riemann-Hilbert Problem 3.2 for Q{z): Find a matrix-valued Junction Q such that: 

1. Q is holomorphic for z G C\rQ. 

2. Asz^ 00, Q{z) =I + (l/z). 

3. Q takes continuous boundary values Q+ and Q- for z & Tq which satisfy the jump relation 
Q+ = Q-Vq where 



(47) 



Vq 



(1 
(1 



koP)^'^^ 



z G (-00,^1) 



and 



(48) 



(0) 

Q 



i?o 

Rq ^ R 



Rn 



zGTo 

zgTi 
z€T2 

ZGT3 



(0) 



Q 



rI 

Ria*/a*o)- 



zGT*o 

z e F* 
^eF^ 

z€Tt 



3.3. Model RHP outside the support. The resulting RHP is an asymptotically stable problem 
in the following sense: as e the jumps remaining on the real axis are non-oscillatory and e 
independent, while along the contours in Tq the jump matrices Vq^ converge to identity both for 

large z and for fixed, nonreal, 2; as e — >■ 0. The convergence in e is uniform on any set bounded 
away from the two stationary phase points ^0 and ^1 where the contours return to the real axis. 
Motivated by these comments we consider the following global model problem which we will later 
prove is a uniformly valid approximation to the solution Q of the RHP defined by (47). Let Uq 
and Ui be, for now arbitrary, fixed size neighborhoods of ^0 and ^1 respectively. We construct a 
parametrix of the form: 



(49) 



'Ao{z) zgUo (cf. 3.3.2) 
P{z)={Ai{z) e (cf. 3.3.3) 

0{z) elsewhere (cf. 3.3.1) 



where the outer model O and the local models ^0 and Ai will be introduced below in the indicated 
subsections. 
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3.3.1. The outer model problem . Away from the real axis , the jumps on UF]^ are exponentially 
small perturbations of identity. By replacing the jump matrices in (47) with their point-wise limits, 
we are lead to the following outer model RHP. 

Riemann-Hilbert Problem 3.3 for the outer model 0{z): Find a 2x2 matrix-valued function 
O such that 

1. O is a hounded function, analytic on C\(—oo,^o]- 

2. 0{z) = I + (l/z) as z^ 00. 

3. For z e (— oo,^o)> O takes continuous boundary values satisfying 0+ = 0_Vo, where 



(50) Vo = 



(l+|roP)2'^3 ^€(-(^,6) 
(l+|ro|2)-3 2 6(^1, Co)- 



This RHP appears in the analysis of the long time problem for NLS [DZ94] . As the jumps are 
diagonal and non- vanishing, one can easily check using the Plemelj formulae [Mus92] that that 
solution is given by 

(51) 0{z) = 5{zy' 

where, 

'Mog(l + |ro(s)|2) , 1 /•«Mog(l + |ro(s)n 



(52) S{z) = exp 



J_ /-^^ log(l + |ro(g)n ^^^ ^ 1 /-^ 
2-jri J_^ s- z 2m i_, 



ds 



The jump matrices for O are discontinuous at and which causes 5 to behave singularly at 
these points. The following lemma describes the nature of the singularities. 

Proposition 2. Define 

f" k(s) 1 
(53) x{z,a):=i ds where k{z) := - — \og{l + \ro{z)\'^) , 

and suppose a and b are fixed real numbers with b < a such that k is analytic in a neighborhood of 
each point. Then near these points: 

i. As z ^ a, X has the uniform expansion, 

X{z, a) = iK{z) \og{z - a) + x{z, a), 

where X is a bounded holomorphic function for z in a neighborhood of a. 

ii. In any suitably small, e independent, neighborhood of b, Mb, the boundary values x+ o-nd X- 
naturally extend as analytic functions to all of Mb and x+i^^ o) ~ X-i^j O') = log(l + |ro(-2)P) 
at each z & Mb. 

Proof. Both results follow from the fact that the weight k in the Cauchy integral defining x is 
analytic at the points a. and b. To prove the first part it is sufficient to observe that the difference 
X — iK,{z) \og{z — a) has a vanishing jump for all z near a. To construct the analytic extensions 
of x± needed near z = b the analyticity of k is used to deform the contour of integration so that 
individually the upper and lower boundary values extend analytically to the opposing half-plane. 
That the extensions satisfy the jump relation follows from the Plemelj formulae. □ 
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Figure 6. Local (exact) jump matrices near the stationary point after intro- 
ducing 6 . 



Expressed in terms of X) = exp[x(2;,^o) +x{z,ii)]', it follows immediately form Prop. 2 
that S has the local expansions 



(54) 



S{z) = {z- ^or^'^S^°\z) near z = ^o, 
S±{z) = {z- ^^r^'^S'l^iz) near z = 6, 



where ^q°' is holomorphic in a neighborhood of and 6^^} are each holomorphic and satisfy 

<^i+/'^i- = 1 + T'o^o any sufficiently small neighborhood of ^i. In particular, observe that (54) 
shows that 5 is bounded at each singularity. The singularity manifests as rapid oscillations as z 
approaches or ^i. 

3.3.2. Local model near z = ^o- Near z = we seek our local model Ao{z) in the form: 

(55) Ao{z)^Aoiz)5{zr. 

This has two advantages. First, it simplifies the matching condition, to match the outer model 0{z) 
the new unknown Aq needs to be asymptotically near identity on the boundary where Ao{z) and 
0{z) meet. Second, as was shown in the construction of the outer model, this factorization removes 
the jump along the real axis. Of course, this comes at the cost of modifying the jumps along the 
other contours. The exact jumps after left multiplication of Q by 6'^^ are given in Figure 6. At 
any fixed distance from the jump matrices arc near identity due to the decay of the exponential 
factors e^*^"/*^ along each ray. The point is a stationary phase point of (cf. (31)), and thus 6o 
is locally quadratic. This allows us to introduce the following locally analytic and invertible change 
of variables 

(56) ^-eo:=-jeo{z)-0o{Co)), or Co = '^^^(^ - 6)- 

The domain Uq on which the local model Aq is defined is selected as follows. We take Uq to be 
any suitably small fixed size neighborhood of such that Uq is bounded away from the contour Fi 

and chosen so that Tinder the map ( = (q{z) its image Co (Wo) is a disk centered at the origin in the 
C-plane (which necessarily has a radius ~ e~^^^). Additionally, we use the freedom to deform the 
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contours Fq and Fs such that the images Co(rfc nZYo), fc = 1, 3 are straight Unes leaving the origin 
at angles 7r/4 and 37r/4 respectively. 

The model inside Uq is constructed by approximating the function ro{z) by its value at ^o, and 
using (54) and (56) to approximate 6 by 

(57) 



Finally, wc make the following change of variables which maps to the C-plane and removes many of 
the constant factors from the jumps: 

- ad <T3 



(58) 



^0 



mo) 



«(?o)/2 



*pc(Co(^);ro(^o)). 



The resulting RHP for the new unknown ^'pc(C) a) is one of the canonical model Riemann-Hilbert 
problems. 

Riemann-Hilbert Problem 3.4 for ^'pc (The Parabolic Cylinder RHP): Given a fixed 

constant a, find a matrix ^'pc((^;a) such that 

1. -^pc is analytic for C S C\{C : arg(C) = ±7r/4, ±37r/4}. 

2. AsC^ oo, *pc = / + *(i)/C + O (l/C') . 

3. On the jump contours, the boundary values satisfy ^pc+ = ^PC-V^pc, where 



(59) 



1 0\ 

^^^-2i«giCV2 1 j 

.0 1 ) 

' -i a* /-1iK,„-iC 

,0 1 



arg(C) = 7r/4 
arg(C) = -7r/4 



arg(C) = 37r/4 
arg(C) = -37r/4, 



1 0^ 

and K = — 5^ log(l + aa*). 
Proposition 3. The RHP for ^pc{C; a) above has a unique, uniformly bounded, solution given by 



(60) 



*pc(C;a) = < 



1 
-a 1 



1 -° 
1 



64 



P(C; a)eK^<^3^-i«<T3 

1 



l + \a\ 

1 a 
1 



arg(C)e(0,7r/4) 

arg(C) e (37r/4,7r) 
|arg(C)|e(7r/4,37r/4) 
arg(C) e (-7r,-3^/4) 

arg(C)e(-7r/4,0) 
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where the function P{(,a), built out of the parabolic cylinder functions -D±iK(-); is 

T^A. (e^C) eTD_,^ (^eT^c) 



(61) P(C;a)=< 



CeC+ 



and i/ie constants k, ^12, anc? ^21 o-tc given by the formulae: 



: C e C- 



1 \/27re''^/*e~'^'*/^ 

K = -—log(l + aa*), /3i2 = ■ . , and /32i = k//5i2- 

ZTT ar(— zk) 

To verify that (60) gives the solution of the RHP for 'I'o one nses the properties of the parabolic 
cylinder functions, _D ,„ (•), (cf [AS64]) to show explicitly that the above formulae satisfies the jump 
and normalization conditions. As this is a somewhat standard model problem we refrain from doing 
so here. The details of the solutions derivation can be found in [DZ94] and [DM08] . 

Later we will need estimates of the error introduced by the local model. On the boundary OUq 

(62) 0-\z)Ao{z) = I + o(e-^/^). 

This follows from (56), the large C, expansion of ^pc (cf. RHP 3.4), and the boundcdness of the 
outer solution 0{z) on dUo. The remaining error introduced by the model lies in the approximation 
of the jumps along Fq, Fa and their conjugates. By direct calculation, one shows that the largest 
contribution to this error is introduced by approximation (57) and consequently along each of these 
rays: 

(63) VQVX^=I + o(e'lHoge). 

3.3.3. Local model near z = ^1. At ^1 the first harmonic 9i is stationary and this forces the lens 
opening contours ri,r2,r3, and their conjugates to return to the real axis at this point. Just as 
in the local model at ^0, the locally quadratic structure of the stationary harmonic motivates the 
definition of a locally invertible analytic change of variables; define C = C,i{z)-i through the relation 

(64) \C,1 := \ {er{z) - ^i(a)) = ^(^ - 6)' + 0{{z- £,,f) . 

We choose the set U\ to be a fixed but sufficiently small neighborhood of ^1 bounded away from 
Wo, such that (64) is analytic and invertible inside, and shaped such that the image Ci(Wi) is a disk 
in the ^-plane. Additionally, we use our freedom to choose the rays Fi, r2, and their conjugates 
such that the images C,i{Tk^Ui), k = 1, 2 lie on the rays arg ( = w/A and 37r/4 respectively. Inside 
Ui we seek a local model of the form 

(65) A^{z) = A^^\z)6{zr^ 

so as to remove the jumps along the real axis exactly from the local problem. The effect of (65) 
on the exact local jumps is to replace Vq by S'^^VqS^'^^; the local contours and exact jumps after 
introducing this S factorization are given in Figure 7. 

Comparing Figures 6 and 7, the jumps near ^1 are more complicated than those near ^o- The 
jump matrices on Fi, F2 and their conjugates involve the full reflection coefficient. As such, at ^1 
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Figure 7. The exact local jumps after conjugating Vq by 6{zy^ are 

[5{z)e-''>^(-y'] ' Vq where Vq are the jumps given above. Inside the shrinking 
disk T> (see (66)) we relabel the contours E/j and define the subsets V^, k = 1,2,3 
in C+ and their conjugates as shown. 



the full multi-harmonic expansion returns to the real axis manifesting as the factor 1 — 7'2g4Lj!y/e 
appearing in the jumps in C+ and the conjugate factor appearing in — and we have to deal 
with every harmonic simultaneously. This is markedly different than the local structure near 
where the jump matrices contained only the single harmonic ^o- However, there arc also important 
similarities, the 6 factors in the off-diagonal entries have a power law singularity of the form (z— ^fe)**^ 
and along the rays Fi, F2, and their conjugates the quadratic decay of a stationary harmonic makes 
the jumps asymptotically near identity at any fixed distance from ^1. These similarities lead one 
to believe that the parabolic cylinder model, RHP 3.4, should be involved in the construction of 
the local approximation. The key to the construction of our model at ^1 is the separation of length 
scales between the onset of asymptotic growth/decay of the locally linear harmonics. Ok, k ^ 1, 
and the locally quadratic harmonic, 61, which wc record in the following elementary proposition. 

Proposition 4. For z <ElAi, the following estimates hold: 

i. For \z — = O (e^/^), the functions e^'^^/^ are hounded independent of e. 

ii. For \z — = O (e), the functions g^'^'*-^*"/'^ are hounded independent of t. 

Hi. Forlmz » e {—Irsiz^e) the function e^^^^l^ (e"'^^"'/'^) is small beyond all orders. 

Proposition 4 motivates the introduction of a shrinking disk inside the fixed sized disk Ui. Define 

(66) V = {z : \z - < ^e}. 

Inside V relabel the jump contours by (^fe)i ^ — lj2, 3, ordered from the real axis to the 
right of ^1 around ^1 with positive (negative) orientation. Additionally, define the real contour 
S4 := {z G 2? : z < ^1}, oriented left to right. Label by 2?o the subset of V bounded by Si and 
S^, and label by V^, 1 < fc < 3, the other subsets of V in C+ ordered coimterclockwisc; let 
denote the complex conjugate sectors. Finally, denote by (V^^) the exact values of the jump 
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5'^^VQd~"^ along T,k i^t)' Figure 7. Our procedure will be to build a model from the parabolic 

cylinder functions as was introduced in the local model at ^o- However, the presence of the linearly 
harmonics requires that inside V the jump matrices be preconditioned for approximation by the 
parabolic cylinder RHP. The price we will pay for this construction are near-identity jumps on dV 
and S4. The preconditioning procedure, which we describe below, is summarized in Figure 8. 
For z € Ui\D, the locally linear harmonics, those terms proportional to 

g*eo/« or e^'^"-"/" in the 

local jumps as expressed in Figure 7, are small beyond all orders (Prop. 4), and they may be set to 
zero without introducing any appreciable error. Note that this approximation replaces the jumps 
on Fq and its conjugate by identity. Eliminating the linear harmonics and employing (54) and (64) 
the remaining jumps are approximated as follows. 

We define the locally analytic and nonzero scalar function 



(67) 



hi = 



6'll\z)6'l^{z)e- 



1/2 



then for z G Ui\D we replace the exact jump 5'^^Vq5 "'^ with K^V^^h-^ , where, 

2 e Fi 



(68) 



1 

-r-o(a)Ci(^)"^'''^«^^e'«^^^^'/' 1 

1 -ra(6)Ci(2)''"'^^^e-<i(^)'/2 
1 

1 



1 



2m(£i)g-iCi(2)V2^ 



-r-o(gi) 

.i+ko(«i)r 



,^^(^)-2i«(6)e»Ci(2)V2 



zgV\ 



zgVo. 




Figure 8. The three steps of approximation in the construction of the model 
problem at ^1 . a) The jumps along the are replaced by the approximation (68) in 
the annular region U\p (bold contours). 6) The exact jumps inside T) are all folded 
onto the real axis to the left of ^1 creating jumps on dD and S4 = r>n{2;<^i}. 
c) The annular approximation are unfolded inside 2? from S4. Dashed contours 
denote jumps which are asymptotically small after the folding and unfolding and 
thus ignored in the construction of the parametrix. 
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Inside V, the linear exponentials cannot be replaced by zero without introducing significant 

errors. Instead, wc use the local self consistency of the jump matrices to fold away the exact jumps 
and then replace them with jumps matching the parabolic cylinder model. To that end we make 
the following change of variables folding all of the jumps inside V onto S4: 

(69) A^^ = A^^^F, 
where F is the piecewise analytic function 

'I z€Vo (I zgUi\V 

(2) 

The new unknown A\ has had the jumps removed from S^, k = 1,2,3, and their conjugates but 
gets new jumps on dT> and S4. A simple calculation shows that the resulting jump along S4 is 

(1 + \ro{z)\'^)~'^^ , while along dV the induced jumps are given by A^'j^} This factorization 

is useful because it has pushed the contribution of the locally linear harmonics off of the contours 
Sfc and onto the disk boundary d{V\DQ) where it is exponentially small away from the real axis. 
Along the arcs &Dz and &D^, which approach the real axis, the induced jumps are independent of 
the linear harmonics: the jump on is given by 

The jump on dV^ is similar. Thus, the locally linear harmonics appear only in the jumps on 
5(2?! U X'2) and its conjugate where their contribution is exponentially small. 

With the exact jumps folded away, we now introduce a second factorization which unfolds onto 
each Ffe nX* the approximation (68). Define 

(70) Af^ = ylf )f/, 
where U is the piecewise analytic function 

(-ro(Ci)Ci(^)-^^^^«^)e<^WV2 1) 



U= < 



' '-ro(6)Ci(^)-'"^«^)e<^(-)V2 iMo 1 I ^ei^s 



^ io 1 / ^ 



'1 ro(6)*Ci(^)^'''^^^^e-'^i(^)V2^ 

.0 1 ) ll^^Ci(^)-^-«^)e<^(^)V2 



elsewhere 



By first folding away the exact jumps (69) and then unfolding the parabolic cylinder model jumps 
(70), the new problem has, by construction, the exact parabolic cylinder jumps along (F^; UF^) nP, 
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fc = 1, 2, 3. The sequence of factorizations induces jumps along the new contours dV and S4, which 
are all near identity. As an example, for z e d'D2, 



FU 



1 




-rp r2 -jgo/eX 
1+roro* X 



(71) 



-ro(a)(l + r-oro*)(5ff)- 
/ 



Se 2e 



ad as 



V 

/ + C(^/eloge) 



(l+r-or2) 



-2iK($i) 



-2i^(«l) 



(2-6)- 



(.-*) 



The calculations in the other sectors are similar with the largest contributing error always being 

■ ^ I 

O (\/eloge) coming from the approximation of {z — ^1)"* by ( v e/^i (^i)Ci I 

Our local model is constructed by simply dropping the near identity jumps along T) and S4. This 
leaves only the jumps on Fi, r2 and their conjugates given by (68) which depend on z only through 
Ci(^). Comparing V^^ with (59) we see that the functions are identical up to the substitution of 
the constant — ro(^i) for ro(^o)- Thus, A^"^^ {C,\{z)) should solve RHP 3.4 up to replacement of the 
appropriate constants. To avoid repeating details, we simply state that our local model in U\ is 
completed by taking 



(72) AX\z) = ^pc (Ci(^),-ro(a)) K"' 

where ^ pc is the function defined in Prop. 3 built from the parabolic cylinder functions. 

3.4. The error matrix, E(z). Proof of Theorem 1, part 1.. Here we prove that the parametrix 
P(^) constructed in the previous subsection is a uniformly accurate estimate of the exact solution 
Q{z) to RHP 3.2 which was derived from the original NLS RHP by explicit transformations. We 
prove this by considering the error matrix Ei^z) defined as the ratio; 

E{,z) = Q{z)p-\z). 

Both Q and the parametrix P are piecewise analytic functions taking continuous boundary values 
on the contours that bound their respective domains of analyticity. As such, E{z) also satisfies a 
Riemann-Hilbert problem with jump relation Ej^ = E-Ve, where 

(73) Ve = P- (VQVf^) PZ\ 

Let F^; denote the totality of contours on which E has a nontrivial jump. Then E satisfies the 

following Riemann-Hilbert problem: 

Riemann-Hilbert Problem 3.5 for the error matrix, E(z). Find a matrix E{z) such that 

1. E{z) is analytic for z S C\Te- 

2. As z 00, E{z) = I + 0{l/z) 

3. For z gTe, E takes continuous boundary values satisfying E_^_ = E_Ve, where Ve is defined by 
(73), 
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Figure 9. The collection, Te, of jump contours for the error matrix, E(z). 



We now shift our perspective and think of E not as defined by (3.4), but as the solution of RHP 
3.5 given the jump matrix Ve which we can calculate explicitly from (47) and (49). The following 
lemma will allow us to establish a uniform asymptotic expansion for E{z). 

Lemma 5. For each {x, t) G K c So compact and £ & No the jump matrix Ve defined by (73) 
satisfies, 

(74) \[z\VE-I)hHr.)=o{e'/Hoge), 
for each sufficiently small e and p= 1,2, or oo. 

Proof. For each {x,t) in a compact subset K of <So, the factorizations defining the RHP for Q can 

be successfully carried through and the distance between Ci and Co is bounded below. For any 
allowable {x,t), the parametrix P is uniformly bonded in the plane, so it suffices to show that 
||(^^(VQVjr^ — 7)11 = O (e^/^loge). The collection of contours Te decouples into three categories: 
the portions of Ffe, A; = 1, 2, 3 and its conjugates lying outside the disks Uq and Ui, these we denote 
F^"*; the portions of each F^ inside Uo and Ui\D together with the real segment E4 we denote F^; 
and finally the loop contours dUo, dU\, and dT> denoted by Fg''^. We verify (74) independently on a 
representative contour of each subset, the derivation on the other contours being similar. Consider 
first the unbounded set of contours F^* as represented by Fg"*. In this region the parametrix is 
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given by the outer model 0{z) which is analytic for z off the real axis. Thus, 
\\^'{Ve - /)IU.(rr) < M\\z'{Vq - /)|Up(rr) 



< M\\z'^-roS-^e'^°/'\ 



Lp(r-*) 



O 



(-'•)■ 



The last line follows from the boundedness of ro and 6 and that, along Fg"*, Im^o increases without 
bound from a least value c > O.The error in the exterior region is small beyond all orders and does 
not contribute to (74). 

On the remaining finite length contours proving that IIVq^p ^ — = O (e^/^loge) implies 

the result for each norm and all moments z^. For the contours in take F|" as representative. 
For z G 2?, Vp = Vg so there is no jump on the inside of T> we need only consider the segment 
w\i\a.D.Ui\D. Using Prop. 4, (47), (67), and (68) we have 



Now let M := max |||/ii5-i||ioo(rin), ||/i^^5||i,o»(rir.) j.Then, 

llVgl/p^ - /||l,oo(rin) 

1 n\ , ^ 

,ro(a)e''^/^Cr''" 



I 2iK(fi) 



— 2iK({i) 

V'-K\) -'^o(?i)Ci 



L«>(ri") 



= 0(^/eloge). 

The final estimate above follows from (64) and the locally quadratic behavior of 6i near ^i. 

Finally consider the disk boundaries OUq, dUi, and 92?. The exact unknown Q is analytic here so 
Vq = I. On the positively oriented boundaries of the two fixed sized disk Vp = h'^^'^ pc{C,k{z); {—^)''^o{ik))h^' 
fc = 0, 1. The e~^/^ scaling in (56) and (64) and the large asymptotics of ^'pc then give 

\\Vf' - Ih^ieu,) = {/ + 0(l/a)} V <0{V^). 

The last contour to consider is the shrinking disk dT>. Using (67)-(70) we have \\Vp^ — I\\L°°(dv) = 
\\d~'^^ {FU — I)5'^^\\Loo(^gx>) = O (\/^foge) where the last equality was previously verified in (71). □ 



Lemma 5 establishes as a small-norm Riemann-Hilbert problem. As such, its solution is given 



by 



f,is)iVEis) - I) 



ds 



s-z 

where /i(s) is the unique solution of (1 — Cy^)/i = /. Here CveI = C-[f{VE — I)] where C_ denotes 
the Cauchy projection operator. In particular, E{z) has the large z expansion 



(75) 



where, 



E^'\x,t) =0(^/eloge). 



Remark. The small norm theory of RHPs as it pertains to this problem can be found in [BC84] 

and [BDT88]. In particular, we need bounds on the Cauchy operators over the contour Fg. For 
the e independent analytic contours in F^; this is a classical result. More care is needed on the 
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e-dependent shrinking contour dV; the boundedness of the Cauchy operators on shrinking circles 
can be found in [KMM03]. 

The proof of part 1 of our main result follows immediately from (75). By inverting the series of 
explicit transformations m t-^ M i-^ Q i-^ E we arrive at a uniform asymptotic expansion of the 
solution m of RHP 2.2. Using (33), (46), and (51) we have 

tp{x, t) — lim 2iz mi2{z; x,t) — O (\/e log e) 

for each fixed {x,t) e Sq. 

4. Inverse problem inside the support, before breaking 







e ^ 






\x-L)_ 
2t 



Figure 10. The regions of growth and decay of the zero harmonic exp {iOi^{z)/e) 
for X € [0,i) . Shaded/Unshaded regions represent regions of decay /growth. The 
dashed line represents the locus of pole accumulation for m{z). 



4.1. Introducing a new pole removing factorization. In this section we being to consider the 
inverse problem for those x £ [0, L), that is, for x inside the support of the initial data. The analysis 
in this regime is made more difficult by the fact that the poles of m in C"*" now lie in the region 
where e*^°/^ is exponentially large, see Figure 10. The exponential growth of e*^"/^ leads us to 
introduce a different pole removing scheme (see (77) below) for x G [0,L). As we will see, the new 
scheme introduces contours with exponentially large jumps which we control by introducing a so 
called g-function transformation which effectively "preconditions" the problem for steepest-descent 
analysis. In this section we show that for each x € [0, L) there exist an order-one finite time Ti{x) 
such that the problem is controlled by the introduction of a genus-zero g-function, i.e., branched 
on a single interval, for all < t < Ti{x), and that the corresponding solution ■ip{x,t) of NLS with 
initial data (2) is asymptotically described by a slowly modulated plane wave. We consider the 
regime x € [0, L) and t > Ti{x) in Section 5. 

To begin our analysis, we introduce a new pole removing factorization which accounts for the 
exponential growth induced by the lowest order harmonic. Let Fi be a semi-infinite contour in C"*" 
leaving the real axis at a point and oriented toward infinity; denote by fl m the region consisting 
of everything to the right of Fi in C+, see Figure 11. For now all we demand of and Fi is that 
the locus of pole accumulation in C+ is contained within ^Im'- 

(76) {0,iq]cnM. 
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In the course of our analysis additional conditions will arise which will serve to more precisely define 
^1 and Fi. 

Remark. For the inverse problem inside the support wc recycle notation from the analysis outside 
the support. Contours and points associated with opening lenses in the RHPs are given the same 
names as before but new definitions, while the functions defining the various matrix factorizations, 
R,Rq, etc., keep their previous definitions. 




Figure 1 1 . The contours Fi , F* and sets CIm and defining the pole removing 
factorization for x E [0, L). The exact shape of the contours Fi and F^, and thus 
the set r^M) ^^M' the point are for now undefined. All we require is that 
flM U fl*j^ contain the locus of pole accumulation (dashed line) and the branch 
F, U F* of = ^/z^ + q^. 



Define 

mR~^Rt) z € flM 

(77) M = < 



m elsewhere. 



The matrices R and Rq are those defined previously, (37) and (43) respectively. Clearly, the new 
unknown M has no poles; the combination mR^^ is analytic in this was the basis of our pole 
removal scheme for x outside the support-while the extra factor Rq has no poles. Additionally, it 
follows from (44) that the nonzero off'-diagonal entry of RRq^ is at leading order proportional to 
gj^i/e g^j^j ^-[^yg unaffected by the exponential growth of e*^"/^. 

Though the new unknown M has no poles it acquires an extra jump in addition to those on the 
lens boundaries. The matrix i?o (43) involved in the pole removing factorization inherits a branch 
cut from I' = a/z^ + which we arc frc^e to choose; we take the cut to be a finite length contour 
F,^ U F* oriented from —iq to iq. Note that as sets F,y and F* are symmetric, but as contours with 
orientation they are antisymmetric. We defer F,^'s exact definition and assume for now only that 
it lies completely within Qm U f2^. The jumps acquired by M along the branch cut are 

(78) i?o-^i?o+ = J) ^eF., 

(79) Rl_Ro\=(l z€Tt. 
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Here we have defined the scalar function 

(80) w := 



ro+ - ro- 



2v+ 
iq 



which we note may be analytically continued away from Vi, fl F*. The RHP for M then follows 
directly from (77). 

Riemann-Hilbert Problem 4.1 for M: Let Tm = M U Ti U U U T* . Find a matrix M{z) 
such that 



(1) M is analytic in C\rM . 

(2) Asz^oo, M{z) =I + (l/z). 

(3) M assumes continuous boundary values andM- onVu satisfying Mj^ 



M-Vm where, 



(WR 

RqRq 



(81) 



V, 



M 




Rq ^r 



z e (-00,^1) 
z e (^1,00) 



zer: 

zeTi 
zgTI. 



4.2. Introducing the g-function. The pole removing factorization introduces jumps on U F* 
which are ill-suited to semi-classical asymptotics. The off-diagonal entries of the jump Vm along F^ 
and F* are proportional to e*^"/*^ and e~*^°/' respectively, which increase exponentially as e — ^ 0. 
To account for these jumps we introduce a scalar 5-function, by making the change of variables 



(82) 



The new unknown N should solve a RHP of the type we have consider thus far, placing the following 
restrictions on g: 

i. g{z) = O (l/z) as z 00. 

ii. There exist a symmetric contour 7U7* such that g is analytic for z € C\(7U7*) with continuous 
boundary values g+ and g- on 7 U 7* . 

iii. g{z)* = g{z*) for all z e C\(7 U 7*). 

Given such a gf-function, the jump matrices for take the following form: 

^Q-''i9+-g-)/<^ ^ei(s++s-)/e 

# #e2^^/' 



Vn = 



2; e 7 U 7* 



2; ^ 7 U 7* 



where the #'s denote the entries of Vm- The g(-function affects the jumps by modifying the expo- 
nential phases in the jump matrices; away from the branch cuts of g the transformation replaces 
the phases 9k in the expansion of the reflection coefficient with 



(83) 



iPk = Ok- 2g. 
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Our procedure is to now construct g so that the resulting RHP for A'' is asymptotically stable, i.e. g 

should remove the exponentially growing components of Vm so that, after additional lens openings, 
the new jump matrix Vn implied by (82) has simple semi-classical limits with controllable errors. 
To construct g, we decompose 7, and by symmetry 7*, into two interlacing sets of contours: bands 
and gaps. Each band is a maximally connected component of 7 satisfying (84a), and each gap is 
a maximal complimentary open interval satisfying (84b); the union of all bands (gaps) in C"*" we 
label 76 (7^) and label the conjugate bands (gaps) 7^ (7g)- 

(84a) Bands: Im (5+ - 5_) = 0. 

f84b) Gaps- / S^^^^) -g+-g-)>0 for z e 7, 

^ I Im (^feo) - 5+ - 5-) < for z G 7*. 

Note that in (84b) we have to select a phase 6'fc(j) which may change from one interval to the next. 
As usual, in the cases where we can successfully construct the g-funetion, it associates with 

the Riemann-Hilbert problem a rational function TZ{x,t) = ^Jfln'lii^ " '^n) and the associated 
Riemann surface of genus G = N — 1. The eventual outer model problem, and thus the leading 
order asymptotics of the solution of NLS, can be described in terms of the Riemann theta functions 
on the Riemann surface TZ{x, t). 

4.3. Genus zero ansatz for x G [0, L) and small times. Having no a priori way to determine 
the number of bands and gaps needed to construct the ^r-function, we proceed instead by ansatz. We 

begin by considering the simplest nontrivial case, that the g-function is analytic for z G C\(r,yUr,y*) 
and the entire contour consists of a single band interval. In this case, we seek g, as the solutions of 
the following scalar RHPs: 

Riemann-Hilbert Problem 4.2 genus zero gf-function and its derivative: Find a scalar 
function g{z) such that 



1. g analytic for z G C\(r^ U T*). 

2. g{z) = 0{l/z) as z ^ CO 

3. 5+ + ,g_ = 6I0 + r?, z&T^yjVl. 

4. 5 = 0((z'±zg)i/2) ^ i^^^iiy 

analytic function near ± iq. 



1. p analytic for z G C\(rj^ U F*). 

2. p{z) = O (l/z^) as z->- 00. 

3. p+ + p- = e'o, zGT.ur*. 

4. p = 0{{z± ig)-^/^) + locally 

analytic function near ± iq. 



Here 77 is a real constant to be determined. The RHP for the density p follows from formally 
differentiating the RHP for g. We begin with the derivative p-problem because it does not depend 
on the unknown constant 77. By integrating the solution of the RHP for p we define a function g 
which a posteriori solves the RHP for g and justifies the formal differentiation. 

Remark. In the typical construction of a g-function, at each band endpoint, aj, the local behavior 
is take to be 5 = O ((z — a^)^/^). However, for square barrier initial data the asymptotic density 
of eigenvalues in the scattering problem (23) is singular: Pq{z) = O ((z ± iq)~^^'^) unlike the usual 
case where eigenvalues vanish as a square root. The loss of regularity in the density motivates the 

O [z ± iqY^^) behavior at the end points. 

Remark. The single band ansatz is motivated by the observation that the initial data ^o{x) is 
locally analytic for x G [0, L). As such the genus zero Whitham equations (4) are locally well-posed 

and so its reasonable to assume that for sufficiently small times the NLS solution, ilj{x,t), is well 
approximated inside the support by a slowly modulating plane wave. In principle we could take 
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,2 



the band to be a subset of U F* with endpoints a{x,t) and a*{x,t). However, in this case one 
can show that the resulting gap condition Im(^o — 9+ — 9-) > for 2; e Fy\7b cannot be satisfied 

which forces a{x,t) = iq. 

RHP 4.2 for p is easily solved using the Plamolj formula: 

f X 1 / n^ ^ 2te^ + (x - L)z + tg^ 

= ^|2^oW-K^) [tz + {x-L)]y 

Any choice of antidcrivativc will satisfy conditions 1,3, and 4 of RHP 4.2 for the gi-function; in 
order for g to tend to zero at infinity, wc define 

(86) g{z) := J p{X) dX = ^Oo{z) - u{z) [tz + {x - L)] + ]^tq 

where to be concrete the path of integration is taken such that it avoids the contour T^, U F* on 
which p is branched. Following (83), g has the effect of replacing the unmodified phases 9q and 9\ 
with 

(87) ^o = 2v[tz + {x- L)]-tq'^ and ipi = 2v [tz -\- {x + L)] - tq^ , 

4.3.1. Satisfying the hand and gap conditions: contour selection. The construction of our ^-function 
is not complete; we have not yet defined the contour U F* on which it is branched. This contour 
is selected by showing that the appropriate band and gap conditions can be satisfied. For the single 
band ansatz we have two such conditions: 

(88a) Band: Im(5(+ - g.) = Im(<^o-) = 0, z 

(88b) Gap: Im(6'i -g+-g_)= Im((^i) > 0, z eV^ 

We list only the conditions for z G C+ as the reflection symmetry implies the correct condition 
in C~. The inequality in (88b) is enforced so that the jumps introduced by (77) along Fi and its 
conjugate are near identity away from the real axis. It is not a typical gap condition since Fi is not 
the complement of a band, but the inequality is of the gap form and, as usual for gap inequalities, 
its failure can be a mechanism for the introduction of a new band in the g-function. 

The band and gap conditions in (88) amount to understanding the topology of the zero level 
curves of Imt^^, fc = 0, 1, and the corresponding open regions of growth and decay of the associated 
exponentials 6*'^'=/'^. Note that the zero level sets of Imipk are independent of how we choose the 
branch cut of v, so using them to determine the branch is not circular. The following lemma 
characterizes the time evolving topology of the zero level sets of each modified phase. 

Lemma 6. Let L{t; b) denote the zero level curve of the function 

Im ipb := Im {2v [tz + b]- tq^) 

for b e K\{0} and v defined by (22). Additionally, define the breaking time 

Ub) = A. 

Then the zero level L(t; b) evolves as follows: 
i. For t = 0, L(0; 6) = M U [-iq, iq\ . 
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ii. For < t < Tc{b), L{t;b) consists of the real axis and two simple non-intersecting arcs: a 

simple finite arc connecting iq to —iq, and an infinite arc which asymptotically a,pproaches the 
line Rez = —b/2t for z large. The finite and infinite arcs lie completely in the half-plane 
Re 62; < and cross the real axis at points zq and z\ respectively such that \zo\ < \zi\. 
Hi. For t > Tc{h), L{t:b) consist of the real axis and two semi-infinite contours each connecting 
one of the branch points ±iq to complex infinity in the same half-plane. In particular, L{t; b) 
does not contain a finite path connecting iq to —iq 

Proof. For t = 0, Iim^b = Imi/ since b is real and nonzero, the topology of L(0; 6) follows immedi- 
ately. For t > the topology is determined from the following observations: 

• The conjugation symmetry (pl{z) = iph^z*) guarantees that M C L{t;b) and allows us to 
consider only z S C"*". Any branch of L{t; b) in C"*" can terminate only at a critical point of 
ipb on the real line, at the branch point iq, or at 00. 

• The square root singularity of </?6 at iq guarantees that exactly one branch of L{t; b) termi- 
nates at iq. 

• The large z expansion ipt, = tz"^ -\- bz -\- t(f /2 + 0{l/z) shows that L{t;b) has exactly 
one branch terminating at 00 in C+ and the branch asymptotically approaches the line 
Re^; = —b/2t for large z. 

• For t > no branch of L{t; b) can cross the imaginary axis, since 



• L{t; b) cannot contain a closed bounded loop if Im</Jb is harmonic inside the loop, since the 
maximum modulus principal would imply that ipb is then identically zero. 

To complete the proof we need to determine the number of real critical points of <^6. A simple 
calculation shows that (pb has exactly two critical points given by 



For < t < Tc{b) thcise are two real points satisfying the given properties in part ii. above. As t 
increases and surpasses the breaking time Tc{b) the two real critical points zq and zi coalesce and 
spit into complex conjugates. Thus, for t > Tc{b) the branches of L{t; b) in C"*" do not have a real 
terminus. Moreover, for t > Tc{b) the branches of L{t; b) no longer pass through the critical points 
since this would necessarily imply the existence of a closed loop in L{t; b) violating the last buUeted 
conditions above. In each case we now connect the termini of the branches of L{t; b) in C"*" without 
violating the above observations. The results for t > follow immediately. □ 

Clearly, the zero level sets of Im (/?o and Ini ipi are those described above with b = x — L and 
b = x-\- L respectively. To satisfy the band condition (88a) we define U F* to be the branch of 
the zero level of Im(/5o connecting —iq to iq. Lemma 6 guarantees this contour exists for {x,t) € <Si, 
where 




0<y<q 

y>q- 




Si := {{x,t) : <x < L,<t < Ti{x)} 



(89) 
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For {x, t) ^ Si the band condition (88a) cannot be satisfied. The upper hmit Ti{x) we call the first 
breaking time, and in the remainder of this section we will always assume that (x,t) £ Si, we will 
address what happens when one moves beyond this first breaking time in Section 5. 

The last condition to be satisfied is the gap condition (88b). For {x,t) g Si the topology of the 
level curve hmpi = falls always into case ii. of Lemma 6. The regions of growth and decay of the 
exponentials e*"^"/^ and 6*"^^/^ follow from continuation of their large z asymptotic behavior; plots 
of these regions for generic values of [x, t) £ Si are given in Figure 12. As the figure illustrates, the 
finite and infinite complex branches of Imipi = do not intersect so the subset of C"*" such that 
Imtpi > consist of the entire region to the right of the infinite complex branch of the zero level 
set except for the closed set bounded by the finite branch of Im (pi and the branch cut of v along 

U F* . To satisfy (88b) we define to be the real critical point of ipi through which the infinite 
branch of Im tpi — passes: 



(90) 



6 = -- 



L 



4t 




1 - 



8i2g2 



{x + L)2 



and take Fi to be any contour in 
Figure 13. 



terminating at laying completely in the set Im ipi > 0, see 




Figure 12. The modified harmonics e*'^''/'^, fc = 0, 1, each separate the plane into 
regions of growth (white) and decay (grey) . The level set Im ipo = (bold lines) 
separating the regions of growth and decay of e'^'^°^'^ consists of three elements: the 
real line, a finite arc connecting ±iq and an unbounded arc which approaches the 
line Imz = — ^'^21!^^ l^-rge z. We take the branch of = \J z^ + along the 
finite branch of Im (p^ — Q which we label F^ in C+ and F* in respectively. The 
level set Im pi — Q consists of similar arcs with the unbounded arc now approaching 
Imz = — for large z. These three arcs and the branch cut Fj, U F* separate 

the regions of growth and decay of e^^^^'^. Finally, the points ^fc, k — 1,2 defined 
by (92) and (90) are the unique points at which the unbounded arcs of Ircupk — 
cross the real line and from which contours are opened. 

4.4. Removing the remaining oscialltions: iV H> Q. Now that the construction of the one 
band g- function is complete the RHP for N follows immediately from RHP 4.1 and (82). 

Riemann-Hilbert Problem 4.3 for N{z): Find a 2 x 2 matrix N such that 
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1. N is analytic in C\rM- 

2. Asz^oo, N{z} = I + 0{l/z). 

3. N assumes continuous boundary values and iV_ on Tm satisfying = N- Vn where, 



(91) 



(g-i(V3o-+«9^)/e 
/ g-»(¥'o-+«9^)/e 

I 





g-i(V3o++tg^)/e 
g-i(¥'o++tg^)/e 



z G (-00,^1) 
Z G (6,00) 



z&Vl 



The g- function transformation has removed the exponentiahy large jumps along r,^ U F*. How- 
ever, the RHP for N still has rapidly oscillatory jumps along the real axis which must be factored 
and moved onto appropriate complex contours. To open steepest descent lenses off the real axis 
we need to understand the regions of growth and decay of the, now modified, first harmonics 
giipo/e g^j^fj 6**^1/^. The topology of these regions follows directly from Lemma 6. The level sets 
Ivsnpk — Q, fc = 0.1 each have one infinite, asymptotically vertical, branch which cross the real axis 
at points which we label ^0 and ^1 respectively; ^1 was given above by (90), following the proof of 
Lemma 6 a similar formula can be given for ^o- 



(92) 



X — L 
At 




These are stationary phase points of the associated harmonics e^'^''^'^. We denote them in this 
way in analogy with the stationary phase points of the same name for the unmodified phases 
9k, k = 0,1 which were described in Section 3.1.1 for x outside the support of the initial data. 
Just as before, these points separate connected unbounded components of the complex plane in 
which the associated harmonic is everywhere either decaying or growing. These regions are plotted 
numerically in Figure 12. 

The oscillatory jumps on the real axis are essentially identical to those in RHP 3.1 up to the 
conjugation by the g-function factor. To deform the oscillations in the current problem onto contours 
on which they are exponentially small we use same factorizations as before (cf. (41)-(45)). Without 
repeating those details we make the following transformation: 

Q{z)=Ne'o^'^''^/''LQ 

(Rq^ 

R-^ {a/aor 



(93) 



-t 



^0 

Rq 

R{a*/al) 



R 



z G flo 

Z 

zefis 
z efil 

z£n*2 

2: G Og 

elsewhere. 
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The regions Clk and the contours Tk bounding them are defined as follows: we take Fq to be a 
contour in lying completely in the region Innpo > and meeting the real axis at ^o! r2 lies 
completely in the region Im(pi < and meets the real axis at ^i; and finally T3 connects the 
stationary phase points ^1 and ^0 and encloses T^, always staying in the region Im<^o < 0- The 
corresponding sets Ofe, A; = 1, 2, 3 are those enclosed by Tk and the real axis, see Figure 13. 




Figure 13. Schematic diagram of the contours Tq and regions flk used to define 
the transformation N ^ Q. 

To calculate the resulting jumps of Q we need the factorizations (41) , (42) , (45) , and the following 
factorization. For z along F^, the jump of Vn can be factored as follows: 

/e-t \ 

1 w-'^e-i'<'°-\ ( -w-ie^*«'\ /l w)-ie-i'^o+ 
1 ) l^w;e-s*9' y VO 1 

Next, recalling the definition of ro(^), (29), and w, (80), we observe the following identity 
which allows one to express the above factorization in terms of the boundary values of W : 

f p-7ifa-+tq') Q \ 

(95) ^ \ , 

This allows us to open a single lens Fa enclosing the band F^, of the ^-function instead of opening 
separate lenses off the real axis and again off F^ which significantly simplifies the construction of 
the ensuing global parametrix. Let F^ = ULo (^fe U 7^) and Fg = F^ U F^ U F* U (-00, ^0] • Using 
this factorization we arrive at: 
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Riemann-Hilbert Problem 4.4 for Q{z): Find a 2x2 matrix valued function Q of z such that: 

1. Q is analytic for z € C\rQ. 

2. Q{z) = I + (l/z) as z^oo. 

3. Q takes continuous boundary values Q+ and Q- on Tq satisfying the jump relation Q+ = Q-Vq 
where, 



(96) 



Vq 







-VJ e 




z&V^ 



z€T* 



l(i + kolT=^ 



and 



Rq 



(97) 



Z&Tq 

zeVi 



Rq ^ R 



(0) 



R 



z e (-00,^1) 

■Rl 

R'iRq^ 
^(aVa5)-3 

i?o 



zen 

z e r*i 

zen 



4- Q admits at worst square root singularities at z = ±iq satisfying the following bounds 

1 \z-iq\-y^ J' 



(98) 



Q{z) = O 
Qiz) = O 



Remark. The new singularity condition in RHP 4.4 is caused by a singularity in (93) defining the 
transformation N ^ Q. In D.^, which contains z = iq, Q = Nc^^^'^'^^^'^Rq^ and observing (94) 
the (2, l)-entry of Rq^ has a square root singularity at iq. We must admit the possibility that Q 
inherits this singularity in its second column. As we encounter such transformations in the rest of 
this paper wc will record the resulting singularity conditions without comment. 

4.5. Constructing a global Parametrix: one band case. The result of the many factorizations 
defining Q is a RHP which is well conditioned for asymptotic approximation. Away from the points 

^0 and 6 where they return to the real axis the jumps (97) along Pq are all near identity, both for 
large z and fixed z as e i 0. The remaining jumps defining Vq behave simply in the semiclassical 
limit. The jumps on the real axis are e-independent and the twist-like jumps on Fj, U F* depend 
on e through an oscillatory constant which wc can deal with explicitly. We proceed, as before, to 
construct a global parametrix P such that the ratio QP~^ is uniformly near identity. Since the 
jumps along Fg are not uniformly close to identity near the stationary phase points 6 6) we 
seek our global parametrix in the form 



(99) 



P{z) 




z€Uo 

z eUi 

elsewhere 
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where each Uk is a fixed size neighborhood of ^/j, k = 1,2, which we will determined in the con- 
struction of the approximation A}.. 

4.5.1. The Outer Mode, 0{z).l. At any fixed distance from the stationary phase points the jumps 
defined by Vq^ in RHP 4.4 are exponentially small perturbations of identity so we ignore them. 
Doing so, we arrive at the following problem for the outer model. 

Riemann-Hilbert Problem 4.5 for 0{z): Find a 2 x 2 matrix valued function O satisfying the 
following properties: 

1. O is an analytic function for z e C\ro, To = Fi, U F* U {—oo, ^o]- 

2. 0{z) = I + 0{l/z) asz^oo 

3. O takes continuous boundary values 0+ and 0_ on Vq away from the points z = ±iq, £,o, and 
^1 . The boundary values satisfy the jump relation 0+ = 0_ Vq where 



(100) Vo 









.w*-'e-i*''' 




zGT^ 



(l + |ro|2)2-3 ^g(_oo,^i) 

Ui + kol')'^^ zei^u^o). 

4- O admits at worst square root singularities at z = ±iq satisfying the following bounds 

^ rnf ^ \z - iq\-^/^ \ 
^(^)=^U \z-iq\-^/^ J' 
(101) ) 
^ ' ' \z + iq\-'/^ 1 



We construct a solution of this outer model by first introdiicing two scalar functions which 
reduces the problem to the class of constant jumps Riemann-Hilbert problems. The first of these 
functions: 



\'27riJ_, 



s- z 2m . _^ s- z , 



we have already encountered; it was introduced in RHP 3.3 to remove the jumps on the real axis 
from the outer model problem. The second function we need is new, and is the solution of the 
following problem: 

Riemann-Hilbert Problem 4.6 for s{z): Find a scalar function s such that 

• s is analytic in C\ (F^ U F*). 

• s{z) = Soo + (1/2) as z ^ 00, where Soo is a constant, bounded for all e. 

• For ^ e F,^ U F* the boundary values s+ and s_ satisfy the relationship 

'u;(z)(5-2(z)e-'*9'/e zeT^ 
w*(^)-i5-2(^)e-'*9'/^ zeTt. 



(102) s+{z)s-{z) = 



At the endpoints o/F^ U F* ; 

i. s{z) = O {{z — iq)^^^) as z — > iq, 
a. s{z) = O ((2 + ig)"^/^) as z ^ —iq. 
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Let US defer, temporarily, the construction of the solutions to RHP 4.6. If such a function s can 
be found, we make the following transformation 

(103) 0{z) = s^^^O'^^\z)s{zy^S{z)''\ 

The new unknown O^^^ then satisfies the following simple, constant jump, Riemann-Hilbert prob- 
lem: 

Riemann-Hilbert Problem 4.7 for O^^^z): Find a 2x2 matrix valued function 0'^^\z) such 
that 

1. 0(1) is analytic in C\ (F^ U T*). 

2. 0(i)(z) ^I + O {l/z) asz^oo. 

3. Away from its endpoints O^^^ takes continuous boundary values O^^ and OL^^ on U F* satis- 
fying the jump realation 

o«(.)=o(^)(.)(; 

4- O^^^ admits at worst 1/A-root singularities at the endpoints z = ±iq. 

This is a standard model problem which often appears in problems related to integrable PDE 
and random matrix theory. One builds the solution by diagonalizing the jump matrix and solving 
the resulting simple scalar problem. A moments work reveals 

( a(z)+a(z)~^ a(z)-a(z)~^ \ 

a{z)-a{z)-^ a{z)+a{z)-^ I ' 
2i 2 / 

where a is defined by, 

(104) a{z)=('^Y' 

\z + iqj 

and the root is understood to be branched along F,y U F* and normalized to approach unity at oo. 

To complete the description of the outer model problem we need to build the solution of RHP 
4.6 which is the subject of the following proposition. 

Proposition 7. Riemann-Hilbert problem 4-6 is solved by 

(105) s{z) = a(.)e-i*«Veexp ( 'M [ ^^-^] . 
Here a is as defined in (104), and using (53) we define 

2(x(^,a) + x(^,Co)) ^eF*. 

where the logarithms are understood to be principally branched. As z ^ oo, s{z) = Soo -\- O {l/z) 
where the constant Sqo satisfies 



flog 


1' 2{z+iq)\ 


[log 




l^2(z-ig) J 



= exp 



-tq - y + «w(a;,t) 



where, 
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Remark. The function co defined by (106) is real and thus represents a slow correction to the fast 
phase. The integrals can be evaluated in terms of special functions to give 

to{x,t) = [Li2 (ro(eo)') +Li2 (ro(6)')] , 

where Li2(z) is the dilogarithm and the {x,t) dependence is contained in the location of the sta- 
tionary phase points and ^i. 

Proof. The Plamelj formulae imply that the given function (105) satisfies the analyticity and jump 
conditions of RHP 4.6. Using classical results on the behavior of singular integrals along the line 
of integration one can show that the Cauchy integral and its boundary values are bounded except 
at the endpoints where at worst it behaves like an inverse square root [Mus92]. However, this 
behavior is balanced by the presence of the factor i^{z) multiplying the Cauchy integral. Thus, 
(105) is everywhere bounded except at z = ±iq where its local behavior is identical to that of a{z). 
Finally, by expanding (105) for large z we get an asymptotic expansion for s at infinity in powers 
of z~^; the leading order term is constant: 




By using contour deformation arguments and exploiting the odd symmetry of v^z) for {z : \z\ > 
^0, 2; € M} and even symmetry for {z : \z\ < ^O) z & R} one shows that satisfies (106). □ 

Note that 0{z) only depends on e through the phase factors e**'' /^"^ appearing in s{z) and Soo 
and is thus bounded independent of e at each point z G C\{iq, —iq}. Near the endpoints z = ±iq, 
0{z) has square root singularities coming from the product 0^^^z)s{zY'^ as described by (101). 
As is expected, and will be proven in Section 4.6, the outer model contributes the leading order 
behavior of the solution of NLS in the semiclassical limit; expanding 0{z) for large z and taking 
the (1, 2)-entry of the first moment we have 

(107) 2i lim Oi2{z; x,t) = iqs^ = gee*9'+*'^(^'*). 

4.5.2. The local model problems. The approximation of Q{z) by the outer model is not uniform 
near z = ^1 and z = ^0 where the jump contours Tq (cf. (96)) return to the real axis. To 
construct a uniformly accurate parametrix we must introduce local models Ao(z) and Ai{z) defined 
on neighborhoods Uq and Ui of the points ^0 and ^1 respectively which account for the local structure 
of the jump matrices and which asymptotically match the outer solution on the boundaries dUo 
and dUi. 

Comparing (96) to (47), the local structure of Vq near ^0 and ^1 when {x,t) G Si, the set 
in which the single band g-function stabilizes the problem, and the local structure of Vq for x 
outside the initial support are essentially identical. In moving from outside the initial support into 
<Si the introduction of the 5-function has replaced the harmonic phases 9k with their modulated 
counterparts (pk = Ok ^ 2g and the stationary phase points ^0 and ^1 have been suitably redefined 
in terms of the (pk, everything else is identical. Thus, we can essentially construct the new local 
models at ^0 and ^1 , mutatis mutandis, from the old models. Without repeating the details of their 
construction wc define the new local models below and refer the reader to Sections 3.3.2 and 3.3.3 
for the additional details. 
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First, define the locally analytic and invertible change of variables C = Ck{z), where, 



(108) 



A: = 0,1. 



Each introduces a rescaled local coordinate at and we choose suitably small fixed size neigh- 
borhoods Uk of ^fe such that the Qk ai'e analytic inside Uk and the images Q = Cfe(^/s) disks in 
the C-plane. Next, using Prop. 2 and (54) we define the nonzero scaling functions: 

1/2 



(109) 



ho 



hi 



mo) 



iK({i 



With the above definitions in hand we define our local models as follows: 



1/2 



(110) 



Mz) = [s-'^^ 0^'\z) s{zr] Ak{z) 6{zr, k = 0,l, 



where the functions Ak are built from the solutions ^pc of the parabolic cylinder local model, RHP 
3.4: 



(111a) ^o(^) = /io-^*Pc(Co(2), ro{Co))ho' 

(111b) ii(z) 



/;.r*Pc(Ci(2), -ro{^i))h^''' X U-^e'3(z)<^^'/ep-l^-tg(z)a,/e_ 



Remark. In addition to the modified phase functions (p^. and the resulting redefinition of the sta- 
tionary phase points ^fe, the pro- factors right multiplying each Ak in (110) arc a new clement of the 
local models when compared to the local models defined in Sections 3.3.2 and 3.3.3. It consists of 
the factors in the outer model (103) which are locally analytic near and ^i; by including them 
as pre-factors, our estimates of the asymptotic error on the boundaries dU^ are identical to the 
previous estimates for x outside the support of the initial data. 

4.6. Proof of Theorem 1, part two. We now consider the error matrix E{z) defined as the ratio: 

(112) E{z) = Q{z)P-\z). 

As both Q{z) and P{z) are piecewise analytic functions whose components take continuous bound- 
ary values on the boundaries of their respective domains of definition, E satisfies its own Riemann- 
Hilbert problem, which we give below. 

Riemann-Hilbert Problem 4.8 for the error matrix, E(z). Find a 2 x 2 matrix E{z) such 
that 

1. E{z) is bounded and analytic for z S C\Te, where Te is the system of contours depicted in 
Figure 9. 

2. Asz-^oo, E{z)^I + 0{l/z). 

3. E takes continuous boundary values E^{z) and E-{z) for z € Te satisfying the jump relation 
E+ = E_ Ve where 

(113) Ve = P- {VqVp^) PZ\ 

All of the properties listed in the RHP for E{z) follow immediately from (112) except for the 

boundedness of E{z) near the endpoints z — ±iq. To show this, first observe that the jump matrices 
of the parametrix P{z) exactly match those of Q{z) along UF* so that Ve, as defined by (113), 
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are identity along r,^ UF*. At worst E{z) has isolated singularities at the endpoints. However, the 

local growth restrictions (98) and (101) dictate that at worst Ej^iz) = O (|z ± zfjj^-'/^) implying 
that the singularities are removable. Thus E{z) is locally bounded near each endpoint. Globally 
boundedness then follows directly from conditions 2 and 3 defining RHP 4.8. 

The key estimate which allows us to prove the second part of our main result, Theorem 1, is 
recorded in the following lemma which establishes E as the solution of a small-norm Riemann- 
Hilbert problem. 

Lemma 8. For {x,t) G K C Si compact and £ G No the jump matrix Ve defined by (113) satisfies, 
(114) \\/{VE-I)\\L.iT.)=o{e'/Hoge), 
for each sufficiently small e and p= 1,2, or oo. 

Proof. The estimates necessary to establish this lemma are identical to those used to prove Lemma 
5 after introducing the modified phases (87) which replace the 9k upon the introduction of the one 
band g-function (86) which remains valid for any {x,t) € iSi. The reader is referred to Lemma 5 
for the details of the estimates. □ 



The small norm estimate in Lemma 8 admits the representation of E{z) as 

^i{s){VE{s) - I) 



ds 



s — z 



where //(s) is the unique solution of (1 — Cy^)/i = /. Here Cve/ = C-[f{VE — I)] where C_ denotes 
the Cauchy projection operator. In particular, E{z) has the large z expansion 



(115) 



E{z) = 1 + 



E(^\x,t) 



where. 



S(i)(.T,t) =0(v^loge) 



Unfolding the sequence of transformations: m i-^ M i-^ N i-^ Q i-^ E, we have the semi-classical 
asymptotic expansion of m and, using (18), we find the leading order behavior of the solution q{x, t) 
of (1) with initial data (2) for each {x,t) e S^. Using (77), (82), (93), (99), and (107) we find 

i^{x,t) = geK9'*+^"(^'*)) +0 (Viloge) . 



4.7. Correction to the phase, evidence of a singular perturbation . We take a moment 

here to remark on the phase correction uj{x,t) defined by (106) that appears in 

V'a.,(:c,f)=geK3'*+-(-.*)), 

the leading order asymptotic behavior of the solution for (x, t) E Si. Wc can think of this phase in 
several ways. Thinking of the asymptotic solution as a slowly modulating plane wave, u) constitutes 
a slow correction to the singular phase. With respect to the fast evolution this is a constant 
phase which the Whitham theory cannot recover. Though the phase is slow it still affects O (1) 
change in the solution, which the Riemann-Hilbert machinery naturally recovers from the scattering 
data. This is a nice result, but can we understand this phase in a more direct way? Altering our 
perspective, we can treat a; as the phase correction in going from the geometric to physical optics 
approximation of the solution in the WKB theory. This leads to an interesting observation. Writing 
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the solution of (1) in the form il){x,t) = A(a;, f)e''^(^'*'/^ we get the following system of equations 
equivalent to focusing NLS 

Pt + {pSx)x = 0, 



(116) 



Pxx 



where p{x,t) := A{x,t)'^. The asymptotic solution ipasy is a correction to the simple plane wave 
solution p = <f and S = q^t. Inserting a regular expansion ansatz: p = q^ + er] + O (e^) and 
S = q^t + eu + (e^), into (116) and differentiating in t to eliminate r] we are led to 

ivtt + q'^Wxx = 0{e). 

Assuming bounded derivatives in the expansion, lu should solve this rcscalcd Laplace equation, but, 
as defined by (106), oj does not satisfy the Laplace equation. To see this, we observe that due to 
the self-similar dependence of and on {x, t) as defined by (90) and (92) we can write w in the 
form uj{x, t) = F (^) + F {^) - coq where 

'■^(f) log(l + |ro(A)|2) 



dX, 



A(C) 



1 



Wo 



loe 



Letting (± 



KA) 

s(i +ko(A )n 

KA) 

{x± L)/t and seeking a solution of this general form, we get 



1 + Vi^VC 



dA. 



1 

t2 



1 d 



c=c+ 



C=C- 



= or F(C) = Ci 



Since C+ and are independent this is only possible if ^ [(C^ 
C2 arctan (^/g). As uj is clearly not of this form, this suggest that the initial assumption that the 
solution admits regular phase and amplitude expansions with bounded derivatives is false. Prom 
the results of [DM05, ER06] we know that the discontinuities exhibit Gibbs phenomena in the 
small-time limit; these arbitrarily small wavelength oscillations are dispersed but not destroyed by 
the evolution. Further, though we know from the semiclassical Riemann-Hilbert analysis that the 
amplitude correction, i], is small, rj = O (^/eloge), it is quite possible that r] exhibits rapid oscil- 
lations with significant derivatives. All of this suggest that the discontinuities generate interesting 
corrections to the phase and amplitude. In particular the phase correction, cj, contributes at first 
order and is detected by the inverse-scattering machinery. We believe that both the phase and am- 
plitude corrections for discontinuous initial data merit further study. Our focus here is to describe 
the regularization of the discontinuities by the semiclassical evolution, so we set these questions 
aside for now, but plan to address them in future work. 

5. Inverse analysis beyond the first breaking time 

In the previous section we showed that for each x € [0, L) the inverse analysis was controlled 

before the breaking time Ti(x), defined by (89), by the introduction of a g- function which a priori 
we assumed to contain a single band interval. The one band ansatz then implied the system of band 
and gap inequalities (88). In this section we now consider the inverse problem for times t > Ti{x). 
The breaking time Ti{x) is characterized by the failure of the band condition (88a), as implied 
by Lemma 6; the lemma also guarantees that the gap condition (88b) continues to be satisfied for 
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times beyond the the first breaking time. The failure of the band condition for a generic choice of 
X G [0, L) as t increases beyond the first breaking time Ti{x) is depicted in Figure 14. In particular, 




Figure 14. The zero level and sign structure of the modified phase IrmpQ intro- 
duced by the one band (/-function as t increases beyond the first breaking time, 
Ti (x) . Shaded regions represent Im ipo > and white regions Im ipo < 0. 



note how the band condition (88a) fails for t > Ti{x): a new gap opens across the real axis. To 
move the analysis of the inverse problem beyond the first breaking time we return to RHP 4.1 for 
the piecewise holomorphic unknown M. In the course of the analysis we will redefine the stationary 
points contours F^,, and introduce a new function which accounts for the newly opening gap. 

5.1. Introducing a new 5- function, the genus one ansatz. As before we introduce our g- 
function by making the change of variables 

(117) N{z) = M(z)e-*»(")'"^/^ 

as always we demand that g satisfy the generic decay and symmetry properties listed below (82). 
Motivated by the appearance of a new gap as discussed above, we seek a new f;-function supported 
on two symmetric bands 7b and 7^ separated by a central gap contour 7gU7* such that 7bU7£, = Fi.. 
We further suppose that the outer endpoints of the bands are fixed at z = ±iq but allow the inner 
endpoints, labelled a and a*, to evolve. Using this ansatz, we seek g as the solution of the following 
scalar Riemann-Hilbert problem: 

Riemann-Hilbert Problem 5.1 for the genus one (^-function Find a scalar function g{z) 
such that 

1. g is an analytic function of z for z £ C\(Fi, U F*). 

2. g{z) = O (1/z) as z ^ 00. 

3. g takes continuous boundary values and g_ on F^ U F* away from the endpoints which satisfy 
the jump relations: 

(.^a^ ( 9+{z)+9-{z)=0o+'n Zejb^lt 

^ ^ { g+{z)-g.{z)^n ze^.ur, 

where rj and Q are real constants. 
4-. As z p, where p is a endpoint of the bands and gaps, g{z) behaves as follows: 

• g{z) = O ((z — pY^"^) for p = iq or p = —iq. 

• g{z) — O [{z — p)"^/^) for p — a or p = a* . 
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If we can construct such a function we define the modified phases as before 
(119) >fik{z)=ek{z)-2g{z), 

but now in terms of the solution g of RHP 5.1. As discussed in Section 4.2, in order for the g- 

function to lead to a scmiclassically stable limiting problem, wc require that it satisfy a system of 
band and gap relations of the form (84). It follows from transformation (117) and RHP 5.1 that 
the necessary band and gap conditions take the form: 

(120a) Bands: Im {(fo-) = for e 76, 

(120b) Gap 1: Im (<^o) > for ^ e jg, 

(120c) Gap 2: Im (v?i) > for ^ e Pi. 

Remark. The above conditions imply, using the reflection symmetry of the (fk, the corresponding 
relations in the lower half-plane. 

Conditions (120a) and (120b) ensure that the g-function docs not introduce large jumps along 
Fj^ur* while the last condition (120c) ensures that the jumps along the contours FiUF^ introduced 
by the pole removing transformation (77) remain small after introducing g. 

5.1.1. Constructing the genus one g-function. We begin by seeking a g-function of the form 

(121) g{z)= f p{X)d\ 

for an unknown density p where the path of integration is any simple contour which docs not 
intersect UP*. The density p must necessarily be symmetric: p(A*)* = /o(A); and solve the 
problem resulting from formally differentiating RHP 5.1 . 

Riemann-Hilbert Problem 5.2 for the genus two density p{z): Find a scalar function p 

such that: 

1. p is an analytic function of z for z € C\(7(, U 7^). 

2. p{z) = O {l/z^) as z^ 00. 

3. On 'fbU'ylp takes continuous boundary values p+ and p- which satisfy the jump condition 

(122) p+{z) + p-{z) = e'o{z), zejbl^lt 

4- As z ^ p, where p is an endpoint of a band or gap interval, p{z) behaves as follows: 



(123) 



p(z) = 0((z-p)-i/2) forp^±iq, 
p{z) — O {{z — pY/"^) for p ~ a, a 



The normalization condition and path restriction guarantee that g \s & well defined function and 
that g{z) = O {\/ z) a,s z ^ 00. The advantage of this representation of g is that the unknown 
constants 77 and 51 arc expressible as particular integrals of the density p. Suppose p is a solution 
of RHP 5.2, then for any z £ and allowable paths for evaluating g+{z) and g-{z) it follows by 
contour deformation that 



(124) g+{z)-g.{z)= p{\)dX = n, zejgU%, 

where Cb is any positively oriented simple closed contoiir in C"*" which encloses jb- Alternatively, 
the allowable paths could be deformed to the conjugate loop, (enclosing 7^ with negative 
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orientation), averaging gives the more symmetric formula: 

(125) n = ^(^J^ pix)dx + J^^p{x)dxy 

which impUes, using the symmetry of p, that f2 is a real constant. 

Remark. Here and elsewhere we make use of the following elementary fact: if f{z) satisfies the 
symmetry /(A*)* = /(A) and C is an oriented contour such that C* = C then 



/(A)dAj = / /(A)*dA*= / /(A*)*dA= / /(A)dA. 
Similarly, if instead C* = —C then: 



^/(A)dA^ =-£/(A)dA. 



Next we consider the jumps of g along the bands and the corresponding constant rj. We have, 

(126) g+{z)+g_{z)=eo{z) + ri, z G 7^ 
where 

(127) r, = -6o{iq) + 2 r p{X)AX. 

J 00 

Repeating the computation for ^; e we have g+{z) + g-{z) = 6q{z) + fj, where, 

-iq 



fi=-9o{-iq) + 2 I p{X) 



dA. 



Observing symmetries, 77 and f] are complex conjugate: fj = rj* . However, a necessary component 
of RHP 5.1 requires that ry be a real constant. By contour deformation, we can express the difference 
as a single integral over the gap: 



V*-V= [ 9'o{X)-2p{X)dX 



"I9 

For generic a. this integral will be nonzero. The vanishing of this integral constitutes a single real 
condition symmetry implies the integral is imaginary on the moving endpoint a. This condition 
has another important implication: in terms of the modified phase lpq^ the vanishing of the above 
condition becomes 

(128) / 6io(A) - 2p(A) dA = / d9?o = 0. 

Jig Jig 

This implies that the endpoints a and a* lie on the same level set of Imt^O) a necessity in order to 
satisfy the band condition (120a). 

We now turn to directly constructing the density p. Motivated by the endpoint behavior and 
the additive jump condition we introduce 

1/2 / * ^ 1/2 



iq J \ z + iq 
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understood to be branched along the bands 7(, U 7^ and normahzed to approach unity as — )• 00. 
Using S{z) and the Plemelj-Sokhotsky formula, it immediately follows that: 

' 27rz 4^7. 5+(A)A-z' 
which can be evaluated by residues yielding the explicit formula 

(129) p{z) = - S{z) [t{2z + a + a*) + {x- L)] . 

Clearly, p satisfies the analyticity and endpoint conditions of RHP 5.2 for any generic choice of a. 
However, for generic a, p{z) = 0{l/z) as z — >■ 00. The quadratic decay at infinity is necessary 
for g to be analytic at infinity and is essential to deriving formulas (125) and (127) for the jump 
constants. In order ior p — O (l/^^^) for large z, a must satisfy the following moment condition: 



where the first equality follows from evaluating the integral by residues. 



-{a + a*) = 0, 



5.1.2. Determining the motion of the a{x,t): self-similar solution of the Whitham system. The 
moment and gap conditions, (130) and (128) respectively, constitute two real conditions on the 
moving endpoints a and a* for each {x, t) in the genus one region. For each i > these conditions 
are equivalent to the intersection of the zero level sets of the functions 

(131a) FM(a,x,t) := ^ (sa^ + 2aa* + 3a* + ''^ ~ (ct + a*), 

nCt 

(131b) FG{a,x,t):= S'(A) [i(2A + a + a*) + (x - L)] dA. 

J a* 

Observing that these functions arc linear — as opposed to affinc functions of x — L and t, it follows 
that for alH > their zero level sets depend only on the self-similar variable 

(132) .:=-^ 

which encodes all of the explicit (x, t) dependence of each function. Geometrically, /i is an angular 
variable parameterizing characteristic lines in {x, f)-space passing through x = L. The exterior 
region. So, with no (/-function, and the genus zero region. Si, are naturally described in terms of 
/Lt as /i < and p > \/2q respcctivcily. The remainder of this section is devoted to the proof of the 
following proposition, which states roughly that for each p in the remaining interval, p e (0, V^q), 
the moment and gap functions (131) uniquely determine a{x,t) := A{p), such that the limiting 
values a(0) and a{y/2q) "match" the exterior and genus one cases. 

Proposition 9. There exist a function A : {0,y/2q) C+ such that by taking a = A{p), 
where p is the self-similar variable (132), the moment and gap equations Fm{A{p'),x' .t') = 
and Fg{A{p'),x' ,t') = 0, with Fm and Fq given by (131), are satisfied for each (x',t') such that 
t' > and p' e (0, \/2g). Additionally, a = A{p) assumes the following limiting values at the 
boundaries of its definition. 

lim A = iq and lim A = -7= 
;*^o fi^V2q V2 
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Proof. Both Fm and Fq are analytic functions of a and a* so the Jacobian necessary to invoke the 
implicit function theorem to find a and a* as functions of x and t is 



(133) J(a,a*):=det 
A straightforward computation using (131) shows that 



dFM dFM 
da da* 

da da* 



(134) 



dFc 


OFm 


r \ -a* 


da 


da 


y„. n{\) 


dFa 


dFM 


X — a 


da* 


da* 





dA 



dA 



where 

7^(A) := ^/{X - iq){X - a)(A + iq){X - a*), 

is cut along the band intervals 7^ and 7^ and normalized such that 7?. ~ A^ for A — )• 00. From this 
it follows directly that 

/1Qr^ *^ / ^.dFu dFu f" dA ^.(a - a*) dFu dFu „, 

(135) J{a,a ) = (" " 5^ |. ^ = ^^^^^ a^^^^(-)' 

where K{m) is the complete elliptic integral of the first kind with parameter 

(136) - = l-f^- 

\a + iq\^ 

Clearly, m G [0, 1] for a G C+, vanishes only for a € M and achieves unity only for a = iq; thus 
the Jacobian is finite and nonzero for each a e C+\{i(7}. Just like the Jacobian, the gap equation, 
Fg, can be expressed in terms of complete elliptic integrals (cf. chpater 22 of [WW96]). Following 
this proc;c'chire, and employing a Landen transformation [BF54] to simplify the form, the gap and 
moment equations Fm = and Fq = are equivalent to the pair of equations 

i (Sa"^ + 2aa* + 3a*'^) + q^ - uia + a*) = 0, 

(137) 4 V y f-v ; , 

\a — iq\^K{m) — Re [{a — iq){a + iq)] E{m) = 0. 

where K and E arc the complete elliptic integrals of the first and second kind respectively. These 
equations together with (136) implicitly define a as a function of n for each a € C~^\{iq}. Using 
these equations we can express a in the convenient form 

(138) a = q (^^/4A{m) - (1 + mA{m))^ + imA{m)^ , 
where we have defined 

(2 - m)E{m) - 2(1 - m)K{m) 
A{m) — ^— T r . 
m^F[m) 

For each m e [0, 1], (138) is well defined and has the following limiting behavior 



q 'iiq , 2\ 

a = —pz + — — m + C m m — >■ 



+ 



72 

a = iq - 2y/l -m + 0{{l-m) log(l - m)) m^l~. 

At these limiting values, we have FM{q/V2,x,t) = 2tq{q — iij^Pi) so fi = \f2q at m = 0, and 
FG{iq,x,t) = —Aiqtii, so /U = at m = 1. This completes the proof. □ 
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Before moving to other considerations, we will show that the endpoints a and iq of the g- 

function describe a special solution for the Ricmann invariants of a slowly modulating one phase 
wave solution of NLS. The level curves Fm = and Fq = ^ define a as a smooth function of both 
X and t which allows one to define the characteristic speed 

/ . N da I da 
(139) c{a,iq):=-- 



dt / dx 



Implicitly differentiating (134) it follows that 



dx J a* Tl(X)""^ ^ dx 



da ia] = Ja-7^(A)— ^ at 



The speed c depends on x and t only through the endpoints iq and a. Using the same procedure used 
to simplify the integrals in the gap equation in the above proof, we can express the characteristic 
speed as 

/, / . X 1^ »N (a — a*){a — iq)K{m) 

(140) c{a,iq) = --{a + a*)- ^ A V ; 



2^ ' {a-iq)K{m) + {iq-a*)E{m)' 

where m is the elliptic parameter defined by (136). Using (139) and (140) we have 

da , ^da 

which is precisely the Whitham equation for the evolution of one Riemann invariant of a slowly 
modulating one phase wave solution of NLS [AAD+94]. The generic one phase modulation is 
described by two conjugate pairs of Riemann invariants labelled A^; each satisfying an evolution 
equation dt\i+Ci{\)dx\i =0, i = 1, . . . , 4. Motivated by the self-similar solution for a described by 
Prop. 9, if we seek solutions of the Whitham system in the form Xi{x, t) := with n as defined 

by (132), the evolution becomes d^^Xi {ci{\) ~ =0- So each invariant Aj is either constant, just 

as the endpoint z = iqin our construction, or they satisfy Cj(A) = /U. Self-similar solutions of this 

type were studied in [Kam97]; there the author shows that setting Ci{\) = ji with q given by (140) 
one arrives at the implicit system of equations (137) which describe the evolution of the endpoint 
z = a. Thus the evolution Ai = iq, \2 = a = A{^), is preciesely a self-similar solution of the 
one-phase Whitham system. 

5.1.3. Topological structure of the zero level sets. With a selected to satisfy the endpoint equations 
(131), the construction of the (^-function will be complete if we can select the band and gap contours 
to satisfy the conditions set out in (120). Those (x,t) for which these conditions can be satisfied 
comprise the region of validity of the genus one (j-function ansatz. Satisfying the band and gap 
conditions amounts to understanding the topology of the zero level sets of Imiyjo and Im(^i. Our 
procedure to determine the structure of the zero level sets mimics the proof of Lemma 6. Let 
Lk{x,t) := {z : Imipkiz} = 0}, the essential facts needed to determine the topology of these sets 
are contained in the following observations 

• The conjugation symmetry ^'^{z) = (pk{z*) guarantees that M C Lk[x,t) and allows us to 
consider only z G C+. 

• Lk{x,t) cannot contain a closed bounded loop if Imipk is harmonic inside the loop as this 
would imply (pf^ is identically zero inside the loop. 
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• Any branch of Lk{x, t) in C+ can terminate only at a critical point of ipk- the branch points 

iq and a, oo, or any other zero of the corresponding density in C+. 
We begin with Lo{x,t), using (119), (121), and (129), we have the convenient representation 



where p is any fixed point in the plane and the path of integration is still understood to avoid 
U r*. Additionally, we have defined the /x-dependent real value 



Simple estimates using (137) and (138) imply that is always positive. Two facts are immediately 
obvious from this representation. First, for t > 0, Lq depends only on the self-similar parameter fx, 
and second, z = £,o is the only critical point of (po not listed above. Using (124) and (126), (po has 
the following local expansions near each critical point, where c is a nonzero constant which depends 
on the expansion point: 

• As z ^ iq, ^o{z) ~ — + c{z — iqY^"^, so one branch of Lq terminates at iq. 

• As ^ — >■ a, {'Po)±{z) ~ —■q^Q. + c{z — a)^/^, so three branches of Lq emerge from a 
separated by angles of 27r/3. 

• As ^; ^ foiz) ~ foiS.o) + c{z — ^o)^j so one branch of Lq in C+ terminates at ^q. 

• As ^; 00, (po(-z) ~ — ^o(^)) so in addition to the real line, one branch of Lq leaves C+ by 
approaching oo along a trajectory asymptotically approaching 'Rez = 

Since the constant terms in each expansion are real, each critical point lies in Lq, and since no 
trajectory leaving a may form a homoclinic orbit the only possibility is that the three trajectories 
leaving a terminate at iq, oo, and respectively. It only remains to determine the topology of 
those connections. Consider the degenerate case /i = /ic = V^?; here the branch point a is real 
and the two bands degenerately form a single continuous path connecting ±iq. Recalling that jic 
parameterizes the first breaking time Ti{x) separating the genus zero and genus one regions, it is 
not surprising that, comparing (85) and (129), the degenerate genus one phase lpq is identically 
equal to its genus zero counterpart along the caustic fi^- Just beyond the caustic, that is for 
jjb = fjbc — h, ?l <C 1, we can view, for any fixed z bounded away from a(/Uc), the genus one phase 
^{)(z) as a small, continuous, perturbation of the genus zero phase. In the genus zero case, we 
showed in the course of proving Lemma 6 that Im(^o(iy) is bounded away from zero for y in any 
closed subset of (0, g). This is robust to small perturbations, and it follows that no trajectory of 
Lq passes through the open imaginary interval between the real axis and iq for /i sufficiently close 
to jjLc- As the trajectories leaving a cannot intersect each other, this completely determines the 
topology oi LQ{x,t) for /x near Hc'- the trajectories leaving a connect, in counterclockwise order, to 
iq, 00, and S^q, see Figure 15. As /i decreases to zero, LQ{x,f?j deforms continuously and thus the 
topology of the level set is preserved. The level set Lo(0 separates C+ into two connected regions, 
in each region the sign of Im (^q is determined by continuation from its limiting behavior for large z. 
The band condition (120a) and the first gap condition (120b) are satisfied by defining 7^ to be the 
unique branch of the level set Im ipQ = connecting a to iq and taking 7^, to be any simple contour 
connecting a to ^0 such that it lies everywhere to the right of Lo{t), that is, in the region of C+ 
where ImLpo > 0. For convenience we choose 7^ such that it approaches a and ^0 along steepest 
descent paths of ipQ, such a choice is depicted in Figure 16. 




p 



(142) 
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lq^ 



ML 



'IN 



\ a/ 

Y 



Figure 15. The connection problem for the level set InnySQ. The left figure shows 
an allowable "twisted" configuration of the branches (dashed lines) connecting each 
of the critical points in C^, while the right figure shows the numerical computed 
branches for the same choice of a. Each possible configuration is topologically 
equivalent and separates C'^ into two connected components on which Imcpo is 
single signed. 
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Figure 16. Numerical computation of Im(^o in C+ for a generic choice of e 
(0, ^/2q). Solid lines represent the zero level set while shaded and unshaded region 
where Im ip^ is positive and negative respectively. The dashed line represents an 
arbitrary choice of the gap contour 7g which we must choose to lie everywhere in 
the region Im <y9o > 0. 



It remains to show that condition (120c) can be satisfied which turns our attention to i^i and 
its associated imaginary zero level set Li{x, t). By definition, (pi = (^o + '^Li>(z) but it also has the 
convenient integral representation 



ipi{z) = ipi{p) + I pi{X)dX 

(143) 



p 



pi(A) = 4i5(A)(A-eo)+4L^ 
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for any fixed base point p in C; as always the path of integration is understood to avoid T^, U T* 
along which v is branched. The two representations lead to the following simple observations: 

• For each z G C^\{iq} such that lm{(po) > 0, lm{(pi) > ALlmviz) > 0. So Li is bounded 
away from Lq except at z = iq. In particular, no branch of Li connects to a. 

• As z — )■ iq, fi{z) ^ —7] + c(z — iqY^'^, so one branch of Li terminates at iq. 

• As z — >■ oo, (p{z) = 2tz'^+2{x+L)z+0 (1/z), so in addition to the real line, one trajectory of 
Li leaves C"*" by approaching infinity along a curve asymptotic to the line Re z — ~(x+L)/2t 

• Near any zero, pk of pi, (pi ^ (pi{pk) + c{z — PkY- If f>i{Pk) is real, four trajectories of Li 
leave pk separated by angles of 7r/2. 

The last points merits further investigation. Clearly, pi(O^) < and lim>_j._tx) Pi(A) = —oo. 
However, for any fixed Aq < 0, hmt^oPi(Ao) = 4LAo/i^(Ao) > 0. Thus, for each sufficiently small! 
t, pi has at least two real, negative zeros. In fact, these two are the only zeros of pi in the complex 
plane; estimates of the quartic equation underlying pi{A) = show that the quartic always has two 
positive real zeros which lie on the other sheet of the Riemann surface associated to pi . Now fix 
p, € [0, V2q], thus fixing the values of a and ^q. For t sufficiently small such that the real zeros exist 
label them and zi, ordered < zi < 0. As t increases the two zeros monotonically approach 
each other eventually reaching a breaking time T2 at which the two zeros coalesce. For times t > T2 
the zeros become complex and ipi has no real critical points. Note that because we solve for T2 
along lines p = const., T2 = T2{x) > Ti{x) for each x e [0, L) (cf. (89)). This breaking time T2{x) 
is uniquely characterized as follows: 

f Pi = 

(144) T2{x) is the unique time such that < has a simultaneous solution. 

For each t > Ti{x), the connection problem is completely determined by the local structure at each 





iq Im ipi 




/ -•.^ 




/ a X 

/ \ 
/ 1 




/ / 



t < T2(X) 





iq Im (pi 
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/ ' 



t = T2(X) 




iq Im (fi 



t > T2(X) 



a \ 



Figure 17. Numerical calculation of limpi as t increases to and beyond T2{x). 
Solid lines indicate where Imipi = and dashed lines the branch cut of (pi Note 
that once t > T2{x), the region Imipi > (shaded region) no longer reaches the 
negative real axis making it impossible to satisfy (120c). 



critical point of (pi. For Ti{x) < t < T2{x) the level set Li{x,t) in C+ consists of the real line, a 
branch connecting ^1 to infinity and a second branch connecting zi to iq. As t increases beyond the 
second breaking time T2{x) the two branches of Li in C"*" meet on the real line and then move into 
the complex plane so that, for t > T2{x), Li has a single branch in connecting iq to infinity, see 
Figure 17. 

The second breaking time t ~ T2{x) limits the validity of the genus-one ansatz for the g- function. 
To satisfy the final condition (120a) the contour Fi must be chosen so that it meets the real axis 
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at a point p such that the region enclosed between Fi and the real interval (p, oo) completely 
encloses the locus of pole accumulation i(0, g), and must lie everywhere in the region Iraipi > 0. 
For Ti{x) < t < T2{x), the condition is satisified by taking Fi such that it leaves the real axis at 
^1, the leftmost real critical point of (pi lying everywhere in the region Innpi > 0. However, for 
t > 72 (x) the set Im^so > pulls away from the negative real axis, so that given any choice of Fi a 
new gap interval emerges across which condition (120c) fails. For t > T2{x) the situation should be 
rectified by replacing the genus-one ansatz with a genus-two ^-function with a new band of support 
across this newly opening gap. However, our focus in this paper is the scmiclassical rcgularization 
of the square barrier so we leave the higher genus transitions as a line of investigation to be pursued 
later. For now we state the following proposition which summarizes the above discussion. 

Proposition 10. Let 

(145) S2 = {{x,t) : X e (0,L), t G {Ti{x), T2{x))} . 

The genus-one ansatz, consisting of the g-function defined by (121) and (129), the band contour"/}, 
connecting iq and a, and the gap contours and Fi exist such that the system of band and gap 
inequalities (120) are satisfied for each {x,t) € S2- 

5.2. Removing the remaining oscillations: N i-^ Q. The completed definition of the ^-function 

results in the following Riemann-Hilbcrt problem for N{z) defined by (117): 

Riemann-Hilbert Problem 5.3 for N(z). Find a 2 x 2 matrix N{z) such that 

1. N{z) is analytic for z € C\Tm- 

2. N{z) =I + (l/z) as z^oo. 

3. For z G Tm N assumes continuous boundary values, 
N+ = N-Vn where, 

'gisad<T3/e (-^ti?) 



and for z G Tm which satisfy 



(146) 



gigad<T3/c (iJ^li?) 



oig ad (73 /e 



(i?ti?-t) 



e-i{<Po-+ri)/e 
''g—i{'Po-+v)/<^ 







«(vo-+r;)/e 



y;g»(vo++f2)/€ giH/e 



z € (-00,^1) 

z e (Ci.oo) 
zeFi 
zgTI 



^g-in/e 
. 



-we 




The jumps along R, Fi, and F^ follow directly from definition (117), while the jumps along the 
bands and gaps follow from (118) and (119). The important observation to make is that, as the 
constants rj and Cl are real, the introduction of the genus-one gf-function has removed the exponential 
growth from the jumps along F^ U F* introduced by the pole removing factorization, replacing it 
with oscillations in the bands and exponential decay in the gaps. 
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To arrive at a RHP which is asymptotically stable in the semiclassical limit we have to now 
introduce factorizations which deform the remaining oscillatory jiimps onto contours on which they 
decay exponentially to identity. These jumps, which lie on the real axis and the bands 7(, and 7^, 
are of the same character as those in RHP 3.1 and RHP 4.3 considered previously in the exterior 
and genus zero cases, respectively. Without repeating the details (c.f. sections 3.2.1-3.2.2), the 
necessary factorizations are identical to those introduced in the preceding cases. To define our 
transformation we introduce contours Ffc and regions flk, A; = 0, 2, or 3 in C+ as follows. Take Fq 
to be a semi-infinite ray leaving 7^ at a fixed point bounded away from a and ^0 lying everywhere 
in the region tpQ > and oriented toward infinity; r2 is a semi-infinite ray lying everywhere in the 
region Im y>i < oriented toward ^1 where it meets the real axis; Fs is a finite contour consisting of 
two pieces, the first oriented from ^1 to a passing over the band 75, and the second oriented from a 
to ^0 passing under 7b: all of T3 lying in the region lm(po < 0. We note as well that the contour Fq 
originates at a point p along the gap contour jg and naturally splits 7^ into two pieces: 7"p oriented 
from p to iq and 7^°'"'^ oriented from ^0 to p. As usual we denote the conjugate contours as 7^^* 
and -y^own* -jijjg corresponding sets Qk are those enclosed by their respective F^ and the real axis, 
see Figure 18. With these definitions in hand, the following transformation defines a new unknown 
Q{z) whose jumps are either near identity or otherwise well approximated by explicit factors for 
which a parametrix can be constructed. 

Define 



Q{z) = Ne'^'^"^/' Lq 



(147) 



Rq ^ 



R-^ [a/aoY 



R, 



-t 



Rq 

R{a*/aiy 

Rl 

I 



2; e fio 
ze9.2 

z 
z 

Z&9.1 
elsewhere. 



Let F^ = ULo (^fe U F^) and Fq = (-00, Co) U F^ U F* U F^. RHP 5.3, the definition of Q, and the 
factorization formulas: (35), (41), (45), and (95) result in the following RHP for the new unknown. 

Riemann-Hilbert Problem 5.4 for Q{z): Find a 2 x 2 matrix Q{z) such that: 



1. Q is analytic for z <E C\Fq. 

2. Q{z) = 1 + (l/z) as z^ 00. 
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3. Q takes continuous boundary values and Q- on Tq satisfying the jump relation Q+ = Q-Vq 

where Vq = 
(148) 



/ -w;-ie'''/A 
[we-'''/'' ) 

( -u;*e'''/' 

\w*-^e-''^/^ 



^€76* 



z € (-00,^1) 
z& (ei,^o) 







^y;g*(¥'o++n)/e giO/e ^ 
^ g-»n/e \ 







Z € 7^°"'" 



^ e 77* 



7 c ^down* 

^ ig 



and V^°^ was previously given by (48). 
'f. Q admits at worst square root singularities at z = ±iq satisfying the following bounds 



(149) 



Qiz)=Oi ^ 



Q{z) = O 



1 \z-iq\-^/^ 



\z - iq\ 



-1/2 



\z + iq\ 
\z + iq\ 



-1/2 I 
-1/2 I 



, z iq, 
, z —iq. 




Figure 18. Schematic diagram of the contours and regions fife in C+ used to 
define the transformation N 1-^ Q. Their counterparts F]^ and in C~ are defined 
by conjugation symmetry. 



5.3. Constructing a global parametrix in the genus one case. The RHP for Q{z) that 
results from the several factorizations above, while defined on an elaborate set of contours, in now 
well conditioned to semi-classical approximations. The jumps along the lens opening contours Fq 
all decay exponentially to identity both as z ^ 00 and as e — >■ 0+ for each fixed z £ Tq away 
from the stationary phase points ^0 and ^1 and the branch points a and a*. The remaining jumps 
along the real axis and the band and gap contours all have simple asymptotic behaviors in the 
semiclassical limit. Along the real axis the remaining jumps are independent of e and in the bands 
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and gaps, the terms in each jump which are not exponentially vanishing are comprised of products 
of e-indepcndcnt factors with constant, e-dcpcndcnt, complex phases. 

We now begin the construction of a global parametrix P with the goal that the error resulting 
from this approximation E = QP~^ is uniformly near identity. The parametrix is necessarily 
picccwisc constructed due to the non-uniformity of the outer approximation near the branch and 
stationary phase points. We seek P in the form: 

Ai z eUi 

Aa Z GUa 

O z elsewhere 

where Uq, Ui, Ua and Ua* are sufficiently small, fixed sized, neighborhoods of ^ij o:, and a* 
respectively. The local problems are all essentially solved: as in the previous cases, the local models 
near and are described by parabolic cylinder functions while the new models introduced near 
the interior endpoints a and a* are described, as we will show, by the standard Airy parametrix. 
The new outer model in the genus one case is completely different than the outer model constructed 
in the previous section for the genus zero problem. 



(150) 



P = < 



5.3.1. The Genus One Outer Model, 0{z). We arrive at the outer model problem by simply drop- 
ping those terms of Vq which are exponentially vanishing for each fixed z as e — t- O"*". This results 
in replacing all of the jumps along Tq with identity and dropping the off-diagonal entries of Vq 
along the gap contours 7g U 7* . The resulting outer model problem is given by: 

Riemann-Hilbert Problem 5.5 for the genus one outer model, 0{z): Find a 2 x 2 matrix 
valued function 0{z) satisfying the following properties: 

1. O is analytic for z G C\To, Tq := U F* U (-00, ,^o]- 

2. 0{z) = 1 + (l/z) as z -> 00. 

3. O takes continuous boundary values 0+ and O- on Tq away from its endpoints. The boundary 
values satisfy the jump relation 0+ = O-Vq where 



(151) 



Voiz) = <^ 



—w{zy 



-l„iT]/e^ 



w(z)e' 



-irj/e 







w*{zy'^e-''^/^ 

''g-in/e \ 







piQ/e 



(l + |ro(z)p)2-3 

1(1 + ^0(^)1^=* 



J 

-w*{z)e''^'^^ 




z e (-00,^1) 
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0{z) 


= o 


(152) 


0{z) 


= o 




0{z) 


= o 



i- O is everywhere bounded except near the endpoints z = ±iq, a, and a* where it satisfies the 
following bounds 

\z + iq\-y^ 1 j' zq, 

\z-p\-^l^ \z-p\-^l^ \ 

\z-p\--l- \z-p\-^l' J' -^V.V = <^ora. 

Remark. The replacement of the jumps in Tq with identity results in the outer solution, 0{z), 
being necessarily singular near each of the endpoints a, a*, and (in addition to the given 
singularities at ±iq). The choice of quarter root singularities at the gap endpoints and bounded 

singularities at the stationary phase points is a choice we make to match the local Airy and parabolic 
cylinder local models which we will introduce in neighborhoods of each pair of points respectively. 

The solution of this problem is given in terms of elliptic theta functions associated with the 
Riemann surface naturally associated to the rational function 



(153) n{z) = -^{z - iq){z - a)[z - a*)[z + iq). 

To be concrete, we will always understand 7?. to be branched along 75 U 7^ and normalized such 
that 'R-{z) ~ z^ for large z. However, to arrive at this solution we must first reduce the outer 
model to the canonical constant jump form by introducing two scalar functions which remove the 
z-dcpcndcnce from the jump matrices. The first of these is the function 8{z) defined by (52). This 
function, as in the previous sections, we use to remove the jumps of 0{z) along the real axis. To 
remove the ^-dependence in the bands we introduce the following scalar Riemann Hilbert problem: 

Riemann- Hilbert Problem 5.6 for s{z): Find a scalar function s{z) with the following prop- 
erties: 

1. s is analytic for z G C\(ri, U T*). 

2. s{z) takes continuous boundary values, s+ and s_, on UF* satisfying 



s+{z)s-{z) = iw{z)5{z) 
(154) s+{z)s-{z) = iw*{z)-^6{z)-^e-^'^'<' 

s+{z)/s-{z) = e-'^/^ 



z € 76, 

z e 7s U7;. 



3. As z^ 00, s{z) = e*P(^) (1 + {1/z)), for some linear function p{z) . 

4- s{z) is bounded and nonzero on any compact set not containing the endpoints z = ±iq where 



s{z) =0 (j^z - iq)^^^^ , as z iq, 
s{z) =o(^{z + i?)"^/^) , as z^ -iq. 



This problem is very closely related to RHP 4.6 which was part of the constant jump reduction 
in the genus zero outer model. 
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Proposition 11. The solution of RHP 5.6 is given by s{z) = so{z)si{z) where 



so{z) = a{z) exp 



7^(z) 

2TTi 



(155) 



r„ur* 



TZ+{X)X-z 



Si{z) = exp 
Here a{z) is given by (104) and 



2e 2m€ 



-ifl dA 



\ 



J(A) = 



log(^ 

VK 



iq) 



2(x(2,a) + x(2;,6)) ze-fb 

zejgUj* 

2(x(z,6) + x(2,Co)) z£jb 



where x{z,a) is defined by (53). Moreover, the polynomial p{z) characterizing the essential singu- 
larity of s at infinity is given by 

1 



(156) 
where 



P{z) =Po{z) + -pi{z) 



Po{z) 



1 

2tt 



r^ur* 



= i I (--A-Rea) 



n 



dA 



Proof. That (155) satisfies tiie jump condition (154) follows from the Plemelj formulae and the 
boundary behavior of TZ{z) and a{z) along their respective branch cuts. The asymptotic behavior 
at infinity and the reality of Pq and pi follow from expanding the Cauchy integral defining ,s(z) and 
the symmetries j(A*)* = -j{X) and {r,y U T*)* = (Fj, U F*)"^ (with respect to orientation of the 
path). Finally, classical estimates from the theory of singular integrals guarantee that the Cauchy 
integrals defining s(z) grow at worst as an inverse square root at each endpoint [Mus92], behavior 
exactly balanced by the rational pre-factor TZ{z). Thus s{z) is bounded and nonzero at each finite 
z except the points ±iq where it inherits the singular behavior of a{z). □ 

Using the functions s{z) and 6{z) we seek the outer model in the form 
(157) 0{z) = I3{z) 0''^\z)s{zy'6{zy', 

where p is the scalar function 

1/4 / * \ 1/4 



(158) 



Piz) = 



iq 



a 



a 



z + iq 



cut along the bands 76U7j and normalized to approach unity for large z. The resulting RHP for the 

new unknown O^^^ (z) is analytic away from the bands where it has constant column permutation 
jumps and square root singularities at half of the branch points. 
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Riemann-Hilbert Problem 5.7 for 0^^\z): Find a 2 x 2 matrix valued function 0^^\z) with 
the following properties: 

1. O''^^ is analytic for z G C\(7f, U 7;^). 

2. Asz^ oo, 0(i)(z) = [I + {1/z)] e-*P(^)'^3 where p{z) is the linear function (156). 

3. 0^^\z) assumes continuous boundary values, o'^\z) and o'^\z), on jb and 7^ which satisfy 
the jump relation 

(159) 0^^\z) = 0^y{z)(^^ J), zGjbU^;. 

4- 0^^\z) is bounded on any compact set not containing the points z = iq or z = a* where it admits 

at worst square-root singularities. 

Riemann-Hilbert problems like O^^^ are well known in the literature of integrable systems and 
random matrices [DIZ97, DKM+99b, KMM03, TVZ04]. We provide the details of its construction 
here for completeness. The construction begins by 'lifting' it onto the genus one Riemann surface 
associated with TZ{z) = \J{z — iq){z — a){z — d){z + iq). We denote this surface by S and label its 
two sheets Ei and S2 arbitrarily. For any non-branch point ^ € C we let P''{z), k = 1,2, denote 
its pro-image on the corresponding sheet. The problem is then lifted onto S by seeking 0^^\z) in 
the form 

(160) 0^'^{z)=[v{P'{z)) , v{P^{z))], 

where v{P) is a single vector-valued fmiction defined for P € E. The jump relation (159) implies 
that the new unknown v{P) is holomorphic away from the branch points and the two pre-images of 
infinity: ooi and 002- Near the branch points z = a and z = —iq the boimded behavior of O^^^z) 
implies that any singularity of v{P) is removable, while the square-root singularities of O^^-* (z) at 
z = iq and z = a* , due to the double-ramification of the branch points, become poles of the function 
v{P). These properties together with the singular behavior at each infinity completely specify the 
function tj{P): 

Problem for v(P): Find a vector- valued function w : E ^ satisfying the following properties: 

• -iTis meromorphic on S\{ooi, 002}, and if (vk) denotes the divisor of the component v^, k = 
1,2, over S\{ooi,oo2} then, 

(161) (Vk) + Pa* + P^g > 0. 

• V is essentially singular at each infinity and 

v(P)e'P^'^ = l + 0(l/z), z^ooi, 

(162) . , , 

v{P)e-'P^'^ = l + 0{l/z), z^0Q2. 

Finding the function v{P) necessary to solve the RHP 5.7 is a classical problem in the study of 
Riemann surfaces whose solution can bo constricted from Baker- Akhiczcr functions. These functions 
are represented here in terms of ratios of the 0-functions corresponding to the Riemann surface S. 
The reader is referred to [FK80] for a general review of the theory of Riemann surfaces and our 
construction of the Baker- Akhiczcr functions follows [DubSl]. As a first step towards finding v, we 
introduce a function / with the following properties: 
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• The component functions fi{P) and f2{P) are meromorphic functions on S such that 

ifl) + Pig + Pa'- 002 >0, 
(/2) + Pig + Pa' - OOi > 0. 

• Each component fk is normaHzed such that /fe(oofc) = 1, k = 1,2. 

The existence and uniqueness of such a function can be proved abstractly, but by direct inspection, 
the functions 



(163) fiAP) = \ 



■^^{z{P)-a){z{P) + iq) 



satisfy the above properties. Here z{P) is the projection of P G S onto C and TZy. is the hfting 
of (153) to the Riemann surface such that TZ-e{P) ~ as P — >■ ooi. Clearly, each component 
fk vanishes at the appropriate infinity and has the correct singularities at Pq* and Piq. Each 
component fk necessarily vanishes at precisely one additional (finite) point in S which we label Pfe. 
Thus, the complete divisor of each component is given by 

{fl)^Pl+^2-P^g~Pa', 
^ ' {f2)=P2+^2-P^g-Pa'■ 

We now consider the ratio ( defined by the component-wise product 
(165) v^f-C 

The resulting unknown ( is the aforementioned Baker- Akhiezer function. That is, the components 
of C have the following properties: 

• Each component (k is meromorphic on S\{ooi, 002} admitting at most a single simple pole 

at the point P^. 

• The local behavior of ^ in a neighborhood of each infinity is given by 



(166) 



where p{z) is the linear function given by (156) and Cfe, k = 1,2 are imknown constants. 

To construct C, we introduce several standard devices on the Riemann surface E. First, we fix 
a homology basis. As S is genus one, the basis consists of only two elements, {a, b}. We take the 
a- loop to lie on both sheets, oriented over the first sheet from (76)- to (j^)^ and the 6-loop we 
take completely on the first sheet as a clockwise loop enclosing 7^ without intersecting or enclosing 
7^, see Figure 19. 

Our choice of homology fixes the basis of holomorphic differentials which consists of exactly one 

element, 

where is a normalization constant chosen so that 

v = 27ri 
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Figure 19. Choice of our homology basis {a, b} on S. Sohd and dotted hnes 
indicate the contour Unes on the first or second sheet respectively. 



We will also need the 6-period of our basis differential, let 

(168) H=L = 2.^I^^ Z/^- 

A standard result in the theory of Riemann surfaces is that Re{H) is strictly negative [FK80]. In 
this case we can say more; the symmetries v{P*)* = v{P) and b* ^ b imply that H is always a 
purely real (negative) number. H allows us to define two important objects. The Riemann constant 
K, we take to be 

(169) K = iTr+^H, 
and we take 

(170) e{w;H)^Y^ exp (^n^H - nw 

as the definition of our Q functions. Q is an entire and even function of w which satisfies the 
automorphic relations 

(171) e{w + 2TTi;H) ^e{w;H), and e{w + H; H) = e-i^-'^Qiw; H). 

The zeros of 8 lie on the lattice w ~ K + 2rmri + nH, n, m £ Z. We choose as the base point on 
S the branch point P = iq, and define an Abel map A ; S — > Jac(E) by 

A{p)= r u. 

J iq 

Finally, let r be the abelian differential of the second kind with double poles at each pre-image of 
infinity, t = tq + e^^Ti: 

(172) T,^^ ntiPT^ -dpk{z{P)), fc = 0,l, 

where p(z) = Po{z) + e~^pi{z) is real linear polynomial (156). The constant Cj- (independent of k) 
is chosen so that 

= 0. 
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Note that at each infinity 

(173) T = dp + o(^^^ , P^ooi, and t = -dp + O ], P ^ 002, 

Let T = To + e~^Ti represent the 6-period of r: 



(174) ^^^fj^' k = l,2, 

and observe that the symmetries t{P*)* = t{P) and b* = b imply that T is necessarily real. 

Lemma 12. The differential n satisfies the relations, ti = —Ct and si{z) exp ^— j Tij = I, 
where the path of integration in the exponential is restricted to the first sheet, Si, and does not 
intersect U F* . 

Proof. Prom (155), Si{z) takes the form si{z) = exp(e~^p(z)), where 

niz) f -in dx 

Using this notation, the second half of the lemma is proved if we can show that p{z) — i /.^ Ti = 0. 
To this end, consider the differential dp. Clearly, dp is holomorphic on C\(7U7*); differentiating the 
jump relations for p wc sec that for z €E 7fc U 7jJ , dp^ + dp^ = and for z € 7,, U 7* , dp+ — dp- = 0. 
Together these facts imply that dp extends naturally to a meromorphic differential on the Riemann 
surface S with singularities only at the two infinities. The difference dp — ir, from (173), is then 
a holomorphic differential. Integrating over an a-cycle, we have for any z G 7h, f^dp — iri = 
p+{z*) + P-{z*) — {p+{z) + P-{z)) = 0. As no nontrivial holomorphic differential can have all its 
a-cycles vanish, the difference must vanish identically. It follows immediately, letting z represent 
an arbitrary point in 7^, that T\ = n = —i dp = —i {p+{z) — p- {z)) = — O. □ 

We now combine these many devices to give Kriechever's formula for the Baker- Akhiezer function. 
Define 

(175) aiP) = N.^^^^^^-^^'^^ 




e{A{p) - A{Pk) 

where the path of integration in the exponential factor is the same as the path in the Abel map 
A{P) and the constant Nk is uniquely chosen to satisfy the normalization condition (166). 

Lemma 13. The function C,{P) defined by (175) is a well defined function from S — >■ for every 
e > and solves the Baker- Akhiezer problem defined by (165) with the parameters N/. given by: 

(176) ^^- e(.To''-t^O) -P(^^° + ^^-^)^^°)' 
where 

(177) Yo= lim IpoiziP)) 

Proof. That C(-P) is well defined follows immediately from the automorphic relations (171) for 6(P) 

and the periods of r. Let us now show that ( is the required Baker- Akhiezer function. Clearly, for 
every choice iV^, ^ 0, (^i^ is meromorphic and its single pole is the unique zero of the 0-function in 
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the denominator. By construction this zero is located at Pfc. One needs only to check that Nk can 
by chosen such that the normalization condition (166) is satisfied. Clearly, one must set 



(178) 



1 



e(A(oofe) - A(Pfc) -K-iT) 



lim exp 



p{z{P)) 



e(A(oofc) - A{Pk) - K) 
which defines provided the right hand side is not identically zero. Recall that 

By Abel's Theorem A{{f)) = so A{Pk) = A{iq) + A{a*) - A{ook'). Now ^(ooi) + ^(ooa) = 0, 
and the evaluation of the Abel map at the branching points iq and a* can be exactly evaluated for 
any concrete choice of paths. Choosing to integrate, for convenience, on the first sheet along the + 
sides of the band and gap contours we have A{iq) = and A{a*) = itt + H/2 = K. Inserting these 
values into (178) and using Lemma 12 and the periodicity relations (171) we have 

1 _ e(iTo - ie-^O) 
Nk ~ QW 

Which is always well defined as Tq and are always real and thus avoids the zeros of the 9- 

function. □ 

We now remove the Riemann surface from the picture to write our outer solution 0{z) in terms 
of integrals lying completely on the complex plane. We do so by exploiting the antisymmetry of 
the differentials: both diff'erentials r] = t and r] = v share the property that 

v = -{ v] ■ 



■exp(-iTo+i(-l)'=+iro). 



/ 



To descend the Abel map onto a function from C to C we make the following restrictions. By j4(z) 
we mean an integral from iq to z lying completely on the first sheet such that the path does not 
intersect the band or gap contours U F* . These restriction on the path eliminate the addition of 
o or 6-cyclcs to the value of the Abel map making A well defined on C. 

RecaUing the definitions (157) of O'^H^) and (160) of v{P) we have the following formula for 
the solution of the outer model problem RHP 5.5: 



(179) 



0{z) 



m 

e{iT) 



So (z) exp i-i j To 5{z) 



where T{z) is the matrix 



nz) = 



' I3{z)+I3{z)-'' e{A{z)-A{oo)-iT) 
2 e(A(z)-A(oo)) 

e(A{z)+A{oo)-iT) 
2 e(A(z)+A(oo)) 



I3{z)-I3{z)-^ e{A{z)+A{oo)+iT) \ 
2 e{A{z)+A(oo)) ' 



fl(z)+P(z)- 
2 



&(A(z)-A{aD)+iT) 
B{A{z)~A(oo)) 



That 0{z) as defined by (179) is a solution of the outer model problem follows directly from 
Prop. 11, Lemma 13, and the observation that \{li{z) + {—iy^'^P{z)~^) are the projections of 
l3{z)fj{P^{z)), j, k e {1, 2}, onto C. To control the error problem later it is necessary to know that 
the outer model remains bounded as e — > 0. As a consequence of Lemma 12, formula (179) for 0{z) 
is e-dependent only through the parameter iT = iTo{x,t) — ie~^Q,{x,t) appearing in the arguments 
of the O functions and as n(x,t) is always real valued, and © is 27ri periodic it follows that the 
outer solution is uniformly bounded as e — >■ 0. 
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As one expects the outer solution encodes the leading order asymptotic behavior of the solution 

ip(x, t; e) of the NLS equation (1) for each {x, t) G ^2 where the gcmis-onc ansatz successfully controls 
the Riemann-Hilbert problem. We record here, for out later use, that leading order behavior: 

2i lim zOi2{z;x,t) = 

(180) _ e(o) e{2A{^) + in-ie-^n) 

{1 ^^"^^e(2^(oo)) e(iTo-ie-iO) 



5.3.2. Construction of the local models. We now describe the construction of the local models near 
the band-gap endpoints a, a* and the stationary points and ^i. At each of these points the 
jump matrices Vq are no longer uniformly near identity which in turn implies that the outer model 
0{z) is no longer a uniform approximation of the solution Q{z) to RHP 5.4 in any neighborhood 
of these points. The local models near and are again constructed from the solution ^'pc 
of the parabolic cylinder problem, RHP 3.4. The new feature in the genus one case are the local 
models near the band-gap endpoints a and a*. The construction of the local model near such 
points is nearly canonical. The local error near these points in controlled by installing a local model 
constructed from Airy functions. 

Let us consider first the real stationary phase points and ^i. The jump matrix Vq, given by 
(148), near these points is essentially the same jump matrix appearing in both the genus zero case 
and the quiescent case for x outside the support (cf. (96) and (47)). The local models at each of 
these points are constructed from parabolic cylinder functions; the reader is referred to Sections 
3.3.2 and 3.3.3 for each model's motivation and additional details. 

Define the pair of locally analytic and invertible functions ( = Ck{z)- 



(181) 



2^0 



— 
1^1 — 



e e 



4t5(A)(A-^o)dA 



A 



AtS{X){X-^o)+4:L—-dX 



Each introduces a rescaled local coordinate at and we choose suitably small, fixed size, 
neighborhoods Uk of such that the Cfe are analytic inside Uk and the images ( = Cfe(^fe) are disks 
in the C-plane. Next, using Prop. 2 and (54) we define the nonzero and locally holomorphic scaling 
functions: 



(182) 



hi 



mo) 



mi) 



((5jj°'(z)s+(z))%-*^''(««)/^ 



1/2 



1/2 



With the above definitions in hand we define our local models as follows: 

Mz) = P{z)0^'\z) Ao{z) s{zrS{zr, 
Ai {z) = p{z)0^^^ {z)s{zy^ ii (z) 5{zy' 
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Figure 20. The jump matrices of Vq near the point z = a. 



where the functions are built from the solution ^'pc(C) ci) of the parabolic cylinder local model, 
RHP 3.4: 

(184a) Ao{z) = h^o'^^PciCoiz), ro+(^o))/io 

(184b) Ai{z) = /if *Pc(Ci(-2), -r-o(a))/ir''' x U-\z)e's^'^'''/'F-\z)e-'s^'^'''/'. 

Here, we recall that F{z) and U{z) arc the folding and unfolding factorizations defined by (69) and 
(70) respectively. The definitions (183), (181), and the large C expansion *pc(C) = I + (C"^) 
imply that on the boundaries Uk we have 

(185) 0-\z)Ak{z) = I + O {V~e) . 

Let us now consider the situation near the band-gap endpoints a and a*. The jump matrices, 
up to orientation of the contours, satisfy the symmetry V{z) = V{z*)^ so it is sufficient to consider 
only the problem at z = a. The local jump matrices of Vq near z = a are depicted in Figure 20. 
Here we have used the relation (94) to write the reflection terms in the jump matrices in terms 
of the single locally analytic and nonzero function w{z). We make no approximations of the local 
jump matrices, the only error introduced by the local model will be on the boundary where it meets 
the outer model. To introduce our local model, let Ua be a neighborhood of a, small enough for 
Fj, to disconnects W„ into two components; let Ua+ {Ua- ) denote the component to the left (right) 
of Fy with respect to orientation. Define the pair of functions on Ua- 

fiz)~< , o ^j, and sgupj^) := 

The continuity of (po and jump relations (118) imply that f{z) is analytic in UaX^b- The function 
/(z) motivates the change of variables 

(186) Uz) := (-^/Wj = (-- 5(A)(A-^o)dAj 

which is a locally analytic and invertible and maps Ua to a neighborhood of the origin in the <J-plane. 

We choose the branch such that Ca(76 HZYq) lies along the negative real axis and use our freedom 
to deform F3 and jg such that | arg(^„(F3 (lUa))] = 27r/3 and arg(^„(7g (lUa)) = 0. We seek our 
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outer model in the form 



(187) 



Ao,{z) = K^{z)^m{Uz)) a2 ti;(z)e-i(''+"^snr.(^)) 



-|<73/2 



where Ka{z) is an analytic pre- factor which we will determine later. The matrix 4'Ai(C) then solves 
the following well known Airy Riemann-Hilbert problem in the C-plane. 

Riemann-Hilbert Problem 5.8 for (The Airy RHP) : Find a 2x2 matrix valued function 
^Ai(C) such that 

1. *Ai is analytic for C € C\ {C : | arg(C)| = 0, 27r/3, tt}. 

2. C -> 



3. assumes continuous boundary values and {"ixi)- on each jump contour satisfying 

the relation (*Ai)+(C) = (*Ai)-(C)14A,(C) where 



This problem, like the parabolic cylinder problem, RHP 3.4, is one of the standard local models 
that emerge in the literature of both inverse scattering and random matrix theory, see [DKM+99b] 
or [DKM+99a]. 

Proposition 14. The solution of RHP 5.8 is given by 



(188) 





$(C)e(if'''-*)-3 



arg(C)e(0,f) 




(189) 



*Ai(C) = { 




) 



^(C)e(if'''-*)-3 



arg(C)e(-f ,0) 



where oj := e^'^'/^ and 
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The solution clearly satisfies the jump conditions along the rays | arg(C)| = 27r/3. The normal- 
ization conditions follows immediately from the asymptotic behavior of the Airy functions [AS64, 
p. 448] and to verify the jump conditons on R one needs the Wronskian relations between the vari- 
ous Airy functions [AS64, p. 446] and the fundamental identity Ai{C) + to Ai{ujC) + lj'^ Ai{w'^Q = 0. 
As the solution is standard in the literature, we leave the verification of these facts to the interested 
reader. 

The last step in the construction of the local model inside Ua is to select the analytic pre-factor 
Ka{z) in order to match the outer model along the boundary diia- From the asymptotic behavior 
(188) if we choose 



Ko,{z) = 0{z) L(^)e-K')+f^«g"r.(^))' 



-^3/2 



0-2 



e 12 



-0-3/4 



1 

-1 



then in light of the e scaling in the definition (186) of Caiz) it follows that 
(190) 0-^{z)Ao,{z)^I + 0{e), z e dU^. 

It remains to show that Ka is analytic in Ua ■ From the definition it is clear that Ka is analytic in 
^a\(76U7g), and by explicitly calculating the jumps along these contours using (151) one sees that 
K^t has identity jumps along these contours. It follows that at worst K,-^ has an isolated singularity 
at z = a. However, the growth condition (152) together with the factor Ca{z)~'^^^'^ imply that at 
most Ka grows as an inverse square root and thus the singularity is removable. This completes the 
construction of the local model at z = a. 

Building the model at z = a* is completely analogous to the construction of the model at 
z = a. We summarize the calculations here. We take as our neighborhood of a* the conjugate 
neighborhood of a: Ua* = ilAaf ■ Then the model inside Ua* is given by 

<T3/2 



(191) 
where 



Aa* {Z) = Ka* {z)^UQa* {z)) -W*{z)-^e 



-f (r)+Osgnp, (z)) 



Qa*{z)= 



(po + n + ^^E^viz) 



and 



Ka*{z) = 0{z) -w*{z)-^e 



- J (17+12 sgnp. («)) 



<^3/2 



e 12 
2^ 



Ca*{z) 



2/3 



-CT3/4 



1 
-1 



Essentially repeating the calculations at z = a, this model exactly removes the jumps of Vq inside 
Ua* while on the boundary we get the same error estimate as in (190): 



(192) 



0-'^{z)Aa*{z)=I + 0{e), zedUa*- 



5.4. The Error RHP, proof of Theorem 1 part three. Wc now show that the paramctrix 
defined by eqs. (150), (179), (183), (187) and (191) is a uniformly accurate approximation of the 
solution Q{z) to RHP 5.4. We do this by proving that the error matrix 

(193) E{z) := Q{z)p-\z) 

satisfies a small-norm Riemann-Hilbert problem. We can then prove the existence of and derive 

an asymptotic expansion for E{z). By unravelling the series of explicit transformations this yields 
an asymptotic expansion for m{z) the solution of the inverse scattering problem, RHP 2.2 which 
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in turn gives an asymptotic expansion for the solution q{x, t) of (1) for each (a;, t) G S2 where the 
genus-one g-function successfully stabilizes the inverse analysis. 

As both Q and P are piecewise analytic whose components take continuous boundary values on 
their respective domains of definition, their ratio, E{z) is also piecewise analytic and satisfies its 
own Riemann-Hilbert problem. 




^1 



Figure 21. The jump contour for the error matrix E in the upper half-plane. The 
boundaries dUi, dU2, and dUa are all fixed sized while the inner loop dV scales like 
e^/^. The jumps in the lower half-plane follow from Schwartz reflection symmetry. 

Riemann-Hilbert Problem 5.9 for the error matrix, E{z). Find a 2 x 2 matrix-valued func- 
tion E{z) such that 

1. E is bounded and analytic for z S C\Fe, where Te is the collection of contours depicted in 
Figure 21. 

2. E{z)^ I + 0{\/z) as 00. 

3. E{z) assumes continuous boundary values E^{z) and E-{z) for z € Te satisfying the jump 
relation E+{z) = E-{z)Ve{z) where 

(194) Ve{z) = P-{z)Vq{z)Vp{z)-'P-{z)-K 

and both Vq and Vp are understood to equal identity where Q or P respectively is analytic. 

All of the properties of E{z) above are a straightforward consequence of (193) except the bound- 
edness of E{z) near the endpoints ±iq. Since the parametrix P{z) has jumps exactly matching 
those of Q{z) along 75 U 7^, E{z) has at worst isolated singularities at ±iq. The local growth con- 
ditions (149) and (152) imply that, at worst, Ejk{z) = O [\z± iq\~^^^). Therefore the singularities 
are removable and E is bounded. 

The following lemma establishes that Ve is uniformly near identity and decays sufficiently fast 
at infinity to admit an asymptotic expansion for large z. 

Lemma 15. For any {x,t) & K c S2 compact and ^ e No the jump matrix Ve defined by (194) 

satisfies 

(195) \\z'iVE-I)\\LH^,)=o[e'/'\oge) 
for every sufficiently small e and p = 1,2, or 00. 
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(196) E{z)=I+^^-^^ + ..., where, \E(^\x,t)\= O (e^'^log^ 



Proof. The jump matrix Ve up to the replacement of the phases Ok with their 5- function conjugated 

counterparts (pk (and the resulting deformation of the contours and stationary phase points) arc 
nearly identical to those in Lemma 5. For (x, t) in any compact K G S2 the endpoints a and a* lie 
an order one distance away from ±iq and each other, and the real stationary phase points ^0 and 
^1 exist and are well separated so we can take the four disks Wo, Ui, Ua, and Ua* to have fixed 
radii depending only on K. Furthermore, the genus one (/-function and the phases (pk satisfy the 
necessary inequalities for the estimates in Lemma 5 to go through with only cosmetic changes away 
from the disks Ua and Ua* ■ The Airy models used to define P inside each disks exactly match the 
jumps of Q, so the only new errors are on the disk boundaries. However, (190) and (192) imply 
that II Ve — /||oo = O (e) for z e oUa U dUct' which is subdominant to the O (e^/^ loge) error terms 
produced in Uo and Ui. □ 

The small norm estimate in Lemma 15 allows us to express the solution E{z) of RHP 5.9 explicitly 

271"* Jte s — z 

where /z(s) is the unique solution of (1 — CvE)l-i = I and Cy^/ := C_[/(Ve — J)] here C_ denotes 
the Cauchy projection operator. In particular, the moment bounds in Lemma 15, justify the large 
z expansion 

E(^Hx,t) 
z 

By unravelling the sequence of explicit transformations m^M^N^Q^E, the asymptotic 
bounds on E{z) give us an asymptotic expansion for 'in. Taking z to infinity along the positive 
imaginary axis (the direction does not affect the answer) we have m{z) = E{z)0{z)e^^^''^'^^^'^ and 
it follows from (18) and (180) that 

9(0) e(2^(oo) + iTp - ie-^n) ^.^^y^ 
'e(2A(oo)) e(iTo-ie-iO) 

gives the leading order term and error bound of the solution of (1) with initial data (2) for each 
{x, t) e S2 in the semi-classical limit. 

Acknowledgements 

We wish to thank Peter Miller for his comments which vastly improved the final draft of this 
manuscript. The authors are supported in part by the National Science Foundation under grants 
DMS-0200749, and DMS-0800979. 

References 

[AAD+94] W. F. Ames, R. L. Anderson, V. A. Dorodnitsyn, E. V. Ferapontov, R. K. Gazizov, N. H. Ibragimov, 
and S. R. Svirshchevskii. CRC handbook of Lie group analysis of diffei-cntial equations. Vol. 1. CRC 
Press, Boca Raton, FL, 1994. Symmetries, exact solutions and conservation laws. 

[AS64] Milton Abramowitz and Irene A. Stegun. Handbook of mathematical functions with formulas, graphs, 

and mathematical tables, volume 55 of National Bureau of Standards Applied Mathematics Series. For 
sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, D.C., 1964. 

[BC84] R. Heals and R. R. Coifman. Scattering and inverse scattering for first order systems. Comm. Pure 

Appl. Math., 37(l):39-90, 1984. 

[BDT88] Richard Bcals, Percy Deift, and Carlos Tomei. Direct and inverse scattering on the line, volume 28 of 
Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 1988. 



N . T XX 0(0) e(2A(oo) + iTo - ie-'^n) j,^ / i /a , \ 



THE SEMI-CLASSICAL LIMIT OF FOCUSING NLS FOR A FAMILY OF NON-ANALYTIC INITIAL DATA 71 



[BF54] Paul F. Byrd and Morris D. Friedman. Handbook of elliptic integrals for engineers and physi- 

cists. Die Grundlehren der mathematischen Wissenschaftcn in Einzeldarstellungen mit besonderer 
Beriicksichtigung der Anwendungsgebiete. Bd LXVII. Springer- Vcrlag, Berlin, 1954. 

[BK99] Jared C. Bronski and J. Nathan Kutz. Numerical simulation of the semi-classical limit of the focusing 

nonlinear schrdinger equation. Physics Letters A, 254(6):325 - 336, 1999. 

[BN67] D. J. Benney and A. C. Newell. The propagation of nonlinear wave envelopes. J. Math, and Phys., 

46:133-139, 1967. 

[Bro96] Jared C. Bronski. Scmiclassical eigenvalue distribution of the Zak;harov-Shabat eigenvalue problem. 

Phys. D, 97(4);376-397, 1996. 
[BroOl] J. C. Bronski. Spectral instability of the semiclassical Zakharov-Shabat eigenvalue problem. Phys. D, 

152/153:163-170, 2001. Advances in nonlinear mathematics and science. 
[DIZ97] Percy A. Delft, Alexander R. Its, and Xin Zhou. A Riemann-Hilbert approach to asymptotic problems 

arising in the theory of random matrix models, and also in the theory of integrable statistical mechanics. 

Ann. of Math. (2), 146(l):149-235, 1997. 
[DKM+99a] P. Dcift, T. Kricchcrbaucr, K. T-R McLaughlin, S. Vcnakidcs, and X. Zhou. Strong asymptotics of 

orthogonal polynomials with respect to exponential weights. Comm. Pure Appl. Math., 52(12):1491— 

1552, 1999. 

[DKM+99b] P. Deift, T. Kriecherbauer, K. T.-R. McLaughhn, S. Venakides, and X. Zhou. Uniform asymptotics for 
polynomials orthogonal with respect to varying exponential weights and applications to universality 
questions in random matrix theory. Comm. Pure Appl. Math., 52(11):1335-1425, 1999. 

[DM05] Jeffery C. DiFranco and Kenneth T.-R. McLaughlin. A nonlinear Gibbs-type phenomenon for the de- 
focusing nonlinear Schrodinger equation. IMRP Int. Math. Res. Pap., (8):403-459, 2005. 

[DM08] M. Dicng and K. D. T. -. McLaughlin. Long-time Asymptotics for the NLS equation via dbar methods. 

ArXiv e-prints, May 2008. 

[Dub81] B. A. Dubrovin. Theta-functions and nonlinear equations. Uspekhi Mat. Nauk, 36(2(218)):ll-80, 1981. 

With an appendix by I. M. Krichever. 
[DVZ94] P. Deift, S. Venakides, and X. Zhou. The coUisionless shock region for the long-time behavior of solutions 

of the KdV equation. Comm. Pure Appl. Math., 47(2): 199-206, 1994. 
[DZ93] P. Deift and X. Zhou. A steepest descent method for oscillatory Riemann-Hilbert problems. Asymptotics 

for the MKdV equation. Ann. of Math. (2), 137(2) :295-368, 1993. 
[DZ94] Percy Deift and Xin Zhou. Long-time behavior of the non-focusing nonlinear schrodinger equation: A 

case study. New Series: Lectures in Mathematical Sciences, (5), 1994. 
[ER06] Adrian Espinola-Rocha. Short time asymptotics of the Manakov System. PhD in Applied Mathematics, 

University of Arizona, 2006. 

[FK80] Hershel M. Farkas and Irwin Kra. Riemann surfaces, volume 71 of Graduate Texts in Mathematics. 

Springer- Verlag, New York, 1980. 

[HT73] Akira Hasegawa and Frederick Tappert. Transmission of stationary nonlinear optical pulses in dispersive 

dielectric fibers, i. anomalous dispersion. Applied Physics Letters, 23(3): 142-144, 1973. 

[Kam97] A. M. Kamchatnov. New approach to periodic solutions of integrable equations and nonlineaj: theory 
of modulational instability Phys. Rep., 286(4);199-270, 1997. 

[KMM03] Spyridon Kamvissis, Kenneth D. T.-R. McLaughlin, and Peter D. Miller. Semiclassical soliton ensem- 
bles for the focusing nonlinear Schrodinger equation, volume 154 of Annals of Mathematics Studies. 
Princeton University Press, Princeton, NJ, 2003. 

[KS02] M. Klaus and J. K. Shaw. Purely imaginary eigenvalues of Zakharov-Shabat systems. Phys. Rev. E (3), 

65(3):036607, 5, 2002. 

[KS03] M. Klaus and J. K. Shaw. On the eigenvalues of ZaJcharov-Shabat systems. SIAM J. Math. Anal., 

34(4):759-773 (electronic), 2003. 
[Lax68] Peter D. Lax. Integrals of nonlinear equations of evolution and solitary waves. Comm. Pure Appl. Math., 

21:467-490, 1968. 

[LM07] Gregory D. Lyng and Peter D. Miller. The iV-soliton of the focusing nonlinear Schrodinger equation for 

N large. Comm. Pure Appl. Math., 60(7):951-1026, 2007. 
[MilOl] Peter D. Miller. Some remarks on a WKB method for the nonselfadjoint Zakharov-Shabat eigenvalue 

problem with analytic potentials and fast phase. Phys. D, 152/153:145—162, 2001. Advances in nonlinear 

mathematics and science. 



72 



ROBERT JENKINS AND KENNETH D. T-R MCLAUGHLIN 



[MK98] Peter D. Miller and Spyridon Kamvissis. On the semiclassical limit of the focusing nonlinear Schrodinger 

equation. Phys. Lett. A, 247(l-2):75-86, 1998. 
[Mus92] N. I. Muskhelishvili. Singular integml equations. Dover Publications Inc., New York, 1992. Boundary 

problems of function theory and their application to mathematical physics, Translated from the second 
(1946) Russian edition and with a preface by J. R. M. Radok, Corrected reprint of the 1953 English 
translation. 

[SY74] Junkichi Satsuma and Nobuo Yajima. Initial value problems of one-dimensional self-modulation of 

nonlinear waves in dispersive media. Progr. Theoret. Phys. Suppl. No. 55, pages 284—306, 1974. 
[TVOO] Alexander Tovbis and Stephanos Venakides. The eigenvalue problem for the focusing nonlinear 

Schrodinger equation: new solvable cases. Phys. D, 146(1-4): 150-164, 2000. 
[TVZ04] Alexander Tovbis, Stephanos Venakides, and Xin Zhou. On semiclassical (zero dispersion limit) solutions 

of the focusing nonlinear Schrodinger equation. Comm. Pure Appl. Math., 57(7):877— 985, 2004. 
[TVZ06] Alexander Tovbis, Stephanos Venakides, and Xin Zhou. On the long-time limit of semiclassical (zero 

dispersion limit) solutions of the focusing nonlinear Schrodinger equation: pure radiation case. Comm. 

Pure Appl. Math., 59(10):1379-1432, 2006. 
[Wlii99] G. B. Whitham. Linear and nonlinear waves. Pure and Applied Mathematics (New York). John Wiley 

& Sons Inc., New York, 1999. Reprint of the 1974 original, A Wiley-Interscience Publication. 
[WW96] E. T. Whittaker and G. N. Watson. A course of modem analysis. Cambridge Mathematical Library. 

Cambridge University Press, Cambridge, 1996. An introduction to the general theory of infinite processes 

and of analytic functions; with an account of the principal transcendental functions. Reprint of the fourth 

(1927) edition. 

[Zak72] V. E. Zakharov. Collapse of Langmuir Waves. Soviet Journal of Experimental and Theoretical Physics, 

35:908-+, 1972. 

[ZS71] V. E. Zakharov and A. B. Shabat. Exact theory of two-dimensional self-focusing and one-dimensional 

self-modulation of waves in nonlinear media. Z. Eksper. Teoret. Fiz., 61(1):118— 134, 1971. 

Department of Mathematics, University of Arizona, Tucson, AZ, 85721, USA 

Current address: Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, USA 

E-mail address: rmjenkinaumich.edu 

Department of Mathematics, University of Arizona, Tucson, AZ, 85721, USA 
E-mail address: mclSmath.arizona.edu 



